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ABSTRACT 

We analyze the one-loop effects (strong and electroweak) on the unconven- 
tional top quark decay mode t —>■ h within the MSSM. The results are 
presented in the on-shell renormalization scheme with a physically well moti- 
vated definition of tan [3. The study of this process at the quantum level is 
useful to unravel the potential supersymmetric nature of the charged Higgs 
emerging from that decay. As compared with the standard mode t b, 
the corrections to t — >■ b are large, slowly decoupling and persist at a size- 
able level even for all sparticle masses well above the LEP 200 discovery range. 
As a matter of fact, the potential size of the SUSY effects, which amount to 
corrections of several ten percent, could counterbalance the standard QCD 
corrections and even make them to appear with the "wrong" sign. Therefore, 
if the charged Higgs decay of the top quark is kinematically allowed -a pos- 
sibility which is not excluded by the recent measurements of the branching 
ratio BR{t —>■ b) at the Tevatron - it could be an invaluable laboratory 
to search for "virtual" supersymmetry. While a first significant test of these 
effects could possibly be performed at the upgraded Tevatron, a more precise 
verification would most likely be carried out in future experiments at the LHC. 



1. Introduction 



Recently, the Standard Model (SM) of the strong and electroweak interactions has been 
crowned with the discovery of the penultimate building block of its theoretical structure: 
the top quark, At present the best determination of the top-quark mass at the 

Tevatron reads as follows : 

= 175±6 Gel/. (1) 

While the SM has been a most successful framework to describe the phenomenology of the 
strong and electroweak interactions for the last thirty years, the top quark itself stood, 
at a purely theoretical level -namely, on the grounds of requiring internal consistency, 
such as gauge invariance and renormalizability-as a firm prediction of the SM since the 
very confirmation of the existence of the bottom quark and the measurement of its weak 
isospin quantum numbers]^]. With the finding of the top quark, the matter content of 
the SM has been fully accounted for by experiment. Still, the last building block of the 
SM, viz. the fundamental Higgs scalar, has not been found yet, which means that in spite 
of the great significance of the top quark discovery the theoretical mechanism by which 
all particles acquire their masses in the SM remains experimentally unconfirmed. Thus, it 
is not clear at present whether the SM will remain as the last word in the phenomenology 
of the strong and electroweak interactions around the Fermi's scale or whether it will be 
eventually subsumed within a larger and more fundamental theory. The search for physics 
beyond the SM, therefore, far from been accomplished, must continue with redoubled 
efforts. Fortunately, the peculiar nature of the top quark (in particular its large mass-in 
fact, perhaps the heaviest particle in the SM!- and its characteristic interactions with the 
scalar particles) may help decisively to unearth any vestige of physics beyond the SM. 

We envisage at least four wide avenues of interesting new physics potentially conveyed 
by top quark dynamics and which could offer us the clue to solving the nature of the 
spontaneous symmetry breaking (SSB) mechanism, to wit: i) The "Top Mode" realiza- 
tion(s) of the SSB mechanism, i.e. SSB without fundamental Higgs scalars, but rather 
through the existence of ti condensates ; ii) The extended Technicolour Models; also 
without Higgs particles, and giving rise to residual non-oblique interactions of the top 
quark with the weak gauge bosons ; iii) The non-linear (chiral Lagrangian) realization 
of the SU{2)l X f/(l)y gauge symmetry 0, which may either accommodate or dispense 
with the Higgs scalars. It can also generate additional (i.e. non-standard) non-oblique 
interactions of the top quark with the weak gauge bosons 0]; and iv) The supersymme- 
tric (SUSY) realization of the SM, such as the Minimal Supersymmetric Standard Model 
(MSSM) P], where also a lot of potential new phenomenology spurred by top and Higgs 
physics might be creeping in here and there. Hints of this new phenomenology may show 
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up either in the form of direct or virtual effects from supersymmetric Higgs particles or 
from the "sparticles" themselves (i.e. the i?-odd |^ partners of the SM particles), in par- 
ticular from the top-squark ( "stop" ) which is the SUSY counterpart of the top quark. Due 
to the huge mass of the latter, one expects that the top-stop system is one of the most 
preferential chiral supermultiplets to which the Higgs sector should couple. Therefore, 
top quark dynamics is deemed to be an ideal environment for Higgs phenomenology and 
a most suitable SUSY trigger, if SUSY is there at all. 

In this paper, we shall focus our attention on the fourth large avenue of hypothetical 
physics beyond the SM, namely on the (minimal) SUSY extension of the SM, the MSSM, 
which is at present the most predictive framework for physics beyond the SM and, in 
contradistinction to all other approaches, it has the virtue of being a fully-fledged Quan- 
tum Field Theory. Most important, on the experimental side the global fit analyses to all 
indirect precision data within the MSSM are at least as good as in the SM; in particular. 



the MSSM analysis implies that rrit = 172 ± 5 |T^, a result which is compatible with 
the above mentioned experimental determinations of m^. 

In the MSSM the spectrum of Higgs-like particles and of Yukawa couplings is far and 
away richer than in the SM. In this respect, a crucial fact affecting the results of our work 
is that in such a framework the bottom-quark Yukawa coupling may counterbalance the 
smallness of the bottom mass, mi, ~ 5 GeV, at the expense of a large value of tan (3 - the 
ratio of the vacuum expectation values (VEV's) of the two Higgs doublets (Cf. Section 2) 
- the upshot being that the top-quark and bottom-quark Yukawa couplings (normalized 
with respect to the SU{2) gauge coupling) as they stand in the superpotential, take on 
the form 

\ _ ht _ rrit \ - _ rUb , , 

Af = — — ;= , Ah = — — ;= . I Z 1 

g V2Mwsmf] g V2Mwcos(3 

Thus, depending on the actual value of tan/5, Af, and At can be of the same order of 
magnitude, perhaps even showing up in "inverse hierarchy": Ab > Xt for tan/? > nit /nib. 
Notice that due to the perturbative bound tan /? < 70 one never reaches a situation where 
\t « Afe. In a sense, A^ ^ A^ could be judged as a natural relation in the MSSM; it can 
even be a necessary relation in specific SUSY-GUT models, e.g. those based on t, h and 
r Yukawa coupling unification [|ll|] , at least at the unification point. Furthermore, one 
expects that if such a relation holds, then it is not just the top-stop system, but also the 
bottom-sbottom chiral supermultiplet that could play a momentous role in the quantum 
physics of the top and bottom quarks. Indeed, since the Higgs sector of the MSSM doubles 
that of the SM, and it comes associated with the fermionic SUSY partners -the so-called 
higgsinos-, one expects that in the limit A;, \t there should occur a very stimulating 
dynamics triggered by the presence of a rich variety of potentially large Yukawa-like 
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interactions formed out of the top/stop-bottom/sbottom-Higgs/higgsino fields. 

A particularly brilliant form of this dynamics, on which we shall focus our attention, 
is revealed through the study of the quantum effects on the non-standard top quark 
decay into a charged Higgs: t —>■ b. This decay, which has already deserved some 
attention in the early literature on the subject [T^, |TB[, is not at all excluded by the 
recent measurements (at the Tevatron) of the branching ratio of the standard top quark 
decay, t b, as will be discussed in more detail in Section 2. The quantum effects 

on t if + b, which we wish to compute in the framework of the MSSM at one-loop, 
can be both strong and electroweak like. Of these the conventional strong corrections 
(QCD) mediated by gluons have already been treated in detail [M. Also the subset of 
strong supersymmetric corrections mediated by gluinos, stop and sbottom squarks, i.e. 



the SUSY-QCD corrections, has been discussed in Ref . [pig] . Here, therefore, we will come 
to grips with the remaining part -as a mater of fact, the largest and most difficult part - of 
the MSSM corrections: namely, the multifarious electroweak supersymmetric corrections 
produced by squarks, sleptons, charginos, neutralinos and supersymmetric Higgs bosons, 
which we shall combine with the total strong (QCD +SUSY-QCD) corrections to obtain 
the MSSM correction. 

In the present studyQ, we will closely follow the systematic pathway adopted in our 
previous treatment of the supersymmetric quantum corrections to the canonical decay 



t W'^ b IjT^, 0. However, because of the Higgs particle in the final state, we have to 
incorporate the details of the renormalization of the Higgs sector of the MSSM, which 
substantially alter the analytical counterterm structure of the t b H^-vertex as compared 
to the conventional t b W^^-vertex. This also has dramatical consequences on the quanti- 
tative side, as we shall see. Moreover, we have to include the quantum corrections from 
the supersymmetric Higgs bosons themselves. In spite of their negligible effect on the 
width of the standard decay t —>■ b [l^ - a fact which was not obvious a priori since 
there are potentially large non-oblique Higgs interactions originating from the Yukawa 
couplings (H) - their impact on the alternative decay t b must also be carefully 

evaluated. In certain circumstances, they can be comparable to the virtual effects from 
the genuine supersymmetric particles (sparticles). 

It is worth mentioning that the one-loop analysis of the unconventional top quark 
decay t —>■ b is motivated not only by its obvious interest on its own as a laboratory to 
test the spontaneous symmetry breaking mechanism beyond the SM, but also as a way to 
characterize the SUSY nature of the charged Higgs emerging from that decay. In fact, we 
shall see that important (~ 50%) SUSY radiative corrections can be obtained in certain 
regions of the MSSM parameter space. In one of these regions, tan/3 ~ mt/mf, is large 
"'^A preliminary presentation of these results was given in Ref. [|l6[ . 
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enough such that Af, ~ Xt. Incidentally, we note that these are typical values of tan/? 
characterizing one of the possible regions of the MSSM parameter space highlighted in 
the literature [^] (other regions have also been exploited [^) in order to alleviate some 
formerly claimed "anomalies" observed in the hadronic decay ratios Rf, and Rc of the Z- 
boson into b b and c c, and in general to try to improve the global fit analyses of the MSSM 
to electroweak precision data. However, most of the fuss about the Z-boson "anomalies" 
seems to have declined; at present, Rc is fully understood within the SM, and there only 
remains a marginal 1.8 a discrepancy in Rb with respect to the SM prediction [p^]. It 
should, nevertheless, be clear that the situation with the MSSM is not bad at all. As 
mentioned above, both the SM and the MSSM can comfortably accommodate the present 
high precision data with basically the same level of statistical significance [|10|. 



Thus the MSSM is at present in very good shape and we are perfectly entitled to deal 
with the decay t b from the SUSY point of view. In particular, we wish to quest 

for generic regions of the MSSM parameter space where the decay rate of t — >■ b is 
competitive with the SM decay t b and at the same time to look for regions where 

it is as much sensitive as possible to virtual supersymmetric effects. For definiteness, in 
this paper we take the point of view that the study of the decay t — > b is worthwhile 
provided that its branching ratio is operative at a level BR{t — > b) > 10%, a condition 
which is not excluded by present Tevatron data (see Section 2). From the theoretical 
point of view, this is fully guaranteed provided that tan/? is large enough (> 30). In these 
conditions, the SUSY-QCD one-loop corrections can typically be in the 50% range and 
the electroweak effects induced by the supersymmetric Yukawa couphngs Xt and Af, can 
be rather large (typically 20%) and so all of them are liable to being measured in future 
experiments at the Tevatron and/or at the LHC We shall see that there are scenarios 
where these electroweak supersymmetric corrections could basically constitute the total 
quantum effect from the MSSM, i.e. the net effect after including the complete (standard 
plus supersymmetric) QCD corrections. 

The paper is organized as follows. In Section 2 the lowest order relations concerning 
the Higgs sector and the top quark decay are given. We also discuss the status of the 
charged Higgs decay of the top quark in the light of the recent data from Tevatron, and the 
prospects for its detection. Section 3 discusses the renormalization of the t b if^-vertex in 
the on-shell scheme with a physically well motivated definition of tan j3. In Section 4 we 
present the full analytical formulae for the one- loop corrected partial width r{t b) 
in the MSSM. The numerical analysis and discussion, as well as the conclusions, are 
delivered in Section 5, where we also comment on the feasibility of measuring the computed 
quantum effects in hadron colliders. Finally, we devote three appendices, respectively on 
SUSY Lagrangians, renormalization details and one-loop functions, in order to make the 
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paper as self-contained as possible. 



2. Lowest order relations and determination of 
BR{t h) from experiment 

In this paper we wish to emphasize the possibility that a charged pseudoscalar, H^, 
involved in a possible unconventional decay of the top-quark, t b, be the charged 

Higgs of the MSSM0. A charged Higgs is necessary in the MSSM since Supersymme- 
try requires the existence of at least two Higgs S't/(2)L-doublets with opposite weak- 
hypercharges to give masses to matter and gauge fields (Cf. Appendix A): 

^i=(fi) iY = -l), ^2=(fi) {Y = +l). (3) 

Because of the SUSY constraints, the structure of the Higgs potential of the MSSM 
constructed out of the two doublets takes on the form [p^ : 



V = mj iHil"^ + ml \H2\'^ -ml^ [^Eij Hi Hi + h.c.) 

+ l{9' + 9"){\H,\^-\H,\f + lg^\HlH,\\ (4) 

where ml, m\, are soft SUSY-breaking masses and g, g' are the SU{2) l x f/(l)y gauge 
coupling constants. After SSB, the physical content of the Higgs sector of the MSSM 
consists of one CP-odd ("pseudoscalar") neutral Higgs, A^, two CP-even neutral Higgs 
bosons, h^,H^, and a charged Higgs boson, H"^. Upon due account of the physical 
gauge sector, the masses of the various spinless bosons are determined in terms of just 
three parameters, which can be chosen to be the two vacuum expectation values (VEV's) 
< i?2 >= V2, < H'l >= vi, giving masses to the top and bottom quarks respectively, and 
one physical Higgs mass. However, due to the SSB constraint 



= v\ + vl = 2M^/g^ = 2-3/2^-^1 ^ (174 GeV)^ , (5) 



where is Fermi's constant in /i-decay, in the end only two parameters suffice to com- 
pletely specify the MSSM Higgs masses at the tree-level. Moreover, since we are interested 
in the decay process t —>■ H^ b, it is natural to take Mh± as the physical input mass rather 
than M^o, as frequently done in other contexts. As the second independent parameter, 
one can take the ratio of the two VEV's: tan/3 = V2/V1. Then, in lowest order, we have 



the relations [^| 



1 



Ko,hO = -(M%+MI± ^{M% + M|)2 - AMIM% cos2 2/5 ) , (6) 



^In the MSSM there are several additional, more exotic, 2-body decays of the top quark and also a 
host of 3-body final states worth studying 123] . 
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where M^o < M^o . It is well-known that these formulas become modified at one-loop |25 
In our case, once Mh± is fixed, the other Higgs masses enter the decay rate of t — b 
only through virtual corrections. Moreover, the renormalization of the masses also induces 
a renormalization of the CP-even mixing angle p5[. Although these one- loop effects are 
necessary to guarantee a neutral Higgs spectrum fully compatible with the phenomeno- 
logical bounds near the tan/5 ~ 1 region, one expects that in practice a one-loop shift 
in the masses and couplings just entails a small (two-loop order) correction to the decay 
rate. We have explicitly checked this fact (Section 5). As for the supersymmetric charged 
Higgs, the present LEP 1.5 bound M^o > GOGeV^pH on the CP-odd state translates 
into a lower limit Mfj± > 100 GeV^, which is not significantly modified by the radiative 
corrections (when M^o is taken as an input) |P5|| . 

The charged Higgs can be, as noted above, very sensitive to bottom-quark interactions. 
Specifically, after expressing the two-doublet Higgs fields of the MSSM in terms of the 
corresponding mass-eigenstates, the interaction Lagrangian describing the tbH^-vertex 
reads as follows : 

^Hbt = l^!* b [mt cot pPR + nih tan/? P^] t + h.c. , (7) 

where Vtb is the corresponding Cabibbo-Kobayashi-Maskawa matrix element, and Pl,r = 
(l/2)(l=i=75) are the projection operators on LH and RH fermions. On the phenomenolog- 
ical side, one should not dismiss the possibility that the bottom-quark Yukawa coupling 
could play a momentous role in the physics of the top quark, to the extend of drastically 
changing standard expectations on top-quark observables, particularly on the top-quark 
width. Of course, this is possible because of the potential tan /^-enhancement of that 
Yukawa coupling. 

In the "a-parametrization" , where the input parameters are (a, M\y, Mz, MH,fnf, ...), 
the coupling g on eq.(0) stands for e/sw, where a = ae.mXl'^ = 0) = e^/in and = 
1 — = 1 — M^/Mz- An alternative framework ( "Gi?-parametrization" ) based on 
the set of inputs {Gp, M]y^ Mz, MH,mf, ...) is also useful, especially at higher orders in 
perturbation theory (Cf. Section 3). At the tree-level, the relation between the two 
parametrizations is trivial: 



V2 2M2 



From the Lagrangian (^ , the tree-level width of the unconventional top quark decay into 
a charged Higgs boson in the G^-parametrization reads: 

X [{m^ + ml- Mjj±){m^ cot^ P + ml tan^ (3) + Amoral] , (9) 
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where 

X'l\l,x\y^) ^ ^[l-[x + ym-{x-yf]. (10) 

It is useful to compare eq.(|^) with the tree- level width of the canonical top quark decay 
in the SM: 

T^^\t^W^h) = f-%')^AV^(l,!4,^) 

V87rV2/ nit mi mf 

x[M^(m? + ml) + {m\ - m^f - 2M^] . (11) 

The ratio between the two partial widths becomes more transparent upon neglecting the 
kinematical bottom mass contributions, while retaining all the Yukawa coupling effects: 



r(o)(t^i/+ h) 

r(o)(t ^ b) 



Af2 , 
i TTTT- 



^ tan^ P + cot' (3 



ml 



(12) 



We see from it that if Mh± is not much heavier than Mw, then there are two regimes, 
namely a low and a high tan/? regime, where the decay rate of the unconventional top 
quark decay becomes sizeable as compared to the conventional decay. They can be defined 
approximately as follows: i) Low tan/5 regime: tan/5 < 2, and ii) High tan/3 regime: 
tan (3 > mt/rrib — 35. The critical regime of the decay t b occurs at the intermediate 

value tan/5 = ^Jm^Jmb ~ 6, where the partial width has a pronounced dip. Around this 
value, the canonical decay t —>■ b is dominant over the charged Higgs decay; more 
specifically, for 3 ^ tan /? < 15 the decay rate of the mode t — * b is basically irrelevant 
as compared to the standard mode: BR{t —>■ b) < 10%. Therefore, a detailed study 
of the quantum effects within that interval is of no practical interest. 

Even though the approximate perturbative regime for tan (3 extends over the wide 
range 

0.5 < tan/? < 70, (13) 

we shall emphasize the results obtained in the phenomenologically interesting high tan jS 
region (typically tan/? > 30). As for the low tan/5 range, while BR{t b) can 

also be sizeable it turns out that the corresponding quantum effects are generally much 
smaller than in the high tan/5 case (Cf. Section 5). Still, we find that in the very low 
tan/5 segment 0.5 ^ tan/5 ^ 1 these effects can be of some phenomenological interest and 
we shall also report on them. As mentioned in Section 1, we do have some theoretical 
motivation to contend that at least one of the two regimes i) or ii) may apply. Therefore, 
it is justified to focus our attention on t ^ b, not only as a possibility on its own, 
but also because there may be realistic situations where it could have a non-negligible 
branching ratio. 
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As a matter of fact, and despite naive expectations, the non-SM branching ratio 
BR{t — s> b) is not as severely constrained as apparently dictated by the recent measure- 
ments of the standard branching ratio at the Tevatron, namely, BR{t h) > 70% . 
To assess this fact, notice that the former result strictly applies only under the assump- 
tion that the sole source of top quarks in pp collisions is the standard Drell-Yan pair 
production mechanism qq —>■ ti\2^. Now, the observed cross-section is equal to the 
Drell-Yan production cross-section convoluted over the parton distributions times the 
squared branching ratio. Schematically, 

(Xobs. = J dqdq a{qq^ti) x \BR{t W+ b)\'^ . (14) 

However, in the framework of the MSSM, we rather expect a generalization of this formula 
in the following way: 

dobs. = J dqdq a{qq ^ tt) X \BR{t ^W+b)]"^ 

+ J dqdq a{qq^ gg) X \BR{g t\)\'^ x \BR{t W+ b)\'^ 

+ j dqdqa{qq^bX) ^\BR{ba^tXi)? ^\BR{t^W+b)\^ + (15) 

where g stand for the gluinos, ti for the lightest stop and ba{a = 1,2) for the sbottom 
quarks. One should also include electroweak and QCD radiative corrections to all these 
production cross-sections within the MSSM. For some of these processes calculations 
already exist in the literature showing that one-loop effects can be important [^, 0. 

It should be clear that the observed cross-section on eq.(|l5l) refers not only to the 
standard bW bW events, but to all kind of final states that can simulate them. Thus, 
effectively, we should substitute BR{t b) in that formula by BR{t ^ Xb), and 

then sum the cross-section over X, where X is any state that leads to an observed pattern 
of leptons and jets similar to those resulting from l^-decay. In particular, X = would 
contribute (see below) to the r-lepton signature, if tan (3 is large enough. Similarly, there 
can be direct top quark decays into SUSY particles that could mimic the SM decay of 
the top quark ||23[] . Notwithstanding, even in the absence of the X contributions, eq.(|l5D 
shows that if there are alternative (non-SM) sources of top quarks subsequently decaying 
into the SM final state, b, one cannot rigorously place any stringent upper bound on 
BR{t b) from the present data. The only restriction being an approximate lower 

bound BR{t 6) > 40 — 50% in order to guarantee the purported standard top 

quark events at the Tevatron |jl|] . Thus, from these considerations it is not excluded that 
the non-SM branching ratio of the top quark, BR{t — >"new"), could be as big as the SM 
one, i.e. ~ 50%. 

Notice that at present one cannot exclude eq.(|T^ since the observed form of the 
conventional t — > b final state involves missing energy, as it is also the case for the 
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decays comprising supersymmetric particles. A first step to improve this situation would 
be to compute some of the additional top quark production cross-sections in the MSSM 
under given hypotheses on the SUSY spectrum. Recently, the inclusion of the qq g g 
mechanism followed by the g ^ tti decay has been considered in Ref . pl| , where it is 
claimed that BR{t ii x?) ^ 50%. By the same token, one cannot place any compelling 
restriction on BR{t —>■ b) from the present FNAL data. In particular, if tan jS is large 
and there exists a relatively light chargino with a non- negligible higgsino component, the 
third mechanism suggested on eg. (p!5|) , namely qq baba followed hj ba txi, could 
also be a rather efficient non-SM source of top quarks. Moreover, if 100 GeV ^ Mh± ~ 
150 GeV, then a sizeable portion of the top quarks will decay into a charged Higgs. Thus, 
if either rrit + mi^ ^ rrig < 300 Gel^ and/or rrit + m^- < < 300 Gel^ , so that at 
least one of the alternative SUSY sources of top quark final states contributing to eq.(|l5D 
is available (and m^,m^^ are not too heavy so that the production cross-section is not 
too phase-space suppressed), then one may equally argue that a large branching ratio 
BR{t b) ^ 50% is not incompatible with the present measurement of the top 

quark cross-sections [||] . This could be most likely the case if the frequently advocated 
SUSY decay t — > Xi is kinematically forbidden. Nonetheless, even if it is allowed, it is 
non-enhanced in our preferential large tan (3 region, in contrast to t b. 

Furthermore, it is worth mentioning that the decay mode t —>■ b has a distinctive 
signature which could greatly help in its detection, viz. the fact that at large tan/3 the 
emergent charged Higgs would seldom decay into a pair of quark jets, but rather into a 
r-lepton and associated neutrino. This follows from inspecting the ratio 

T{H+^T+i^r) _ I fmr\^ tan2/3 



T{H+^cs) 3\mJ (m2/m2)tan2/3 + cot2/3 



^ -( — ) > 10 (for tan/3 > JmJms > 2) , (16) 
3 Knis/ * 

where we see that the identification of the charged Higgs decay of the top quark could be a 
matter of measuring a departure from the universality prediction for all lepton channels. 
In practice, r-identification is possible at the Tevatron p2|, p3| ; and the feasibility of 
tagging the excess of events with one isolated r-lepton as compared to events with an 
additional lepton has also been substantiated by studies of the LHC collaborations |3^ . 
The experimental signature for ti b b would differ from ti —>■ 6 6 by an 

excess of final states with two r-leptons and two b-quarks and large missing transverse 
energy. 

A preliminary study in this direction by the CDF collaboration at the Tevatron has 
been able to exclude a large portion of the (tan /?, Mj:^±)-plane characterized by tan P > 60 
and Mh± below a given value which varies with tan p. For extremely high tan /? > C(IOO), 
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the uppermost excluded mass region is Mh± ~ 140 GeV . However, within the interval 
tan/? = 60 — 80, the allowed upper limit on Mh± varies very fast with tan/5. In particular, 
the MSSM permissible values Mh± ~ 100 GeV (compatible with M^o > 60 GeV) are not 
manifestly excluded for tan/3 equal or below the perturbative bound tan/3 = 70, eg. (p!3D . 
We shall nevertheless err on the conservative side and assume that tan /3 < 60 throughout 
our analysis. Thus, as far as the high tan (3 regime is concerned, we will for definiteness 
optimize our results in the safe, and phenomenologically interesting, high tan j3 segment 

30<tan/3<60. (17) 

To round off the r-lepton business, it has recently been shown that it should be fairly 
easy to discriminate between the H^-daughter r's and the if^-daughter r's by just taking 
advantage of the opposite states of r polarization resulting from the and decays; 
the two polarization states can be distinguished by measuring the charged and neutral 
contributions to the 1-prong r-jet energy (even without identifying the individual meson 
states) m, |7|. 



In short, there are good prospects for detecting the decay t — > b, if it is kinemati- 
cally accessible. Unfortunately, on the sole basis of computing tree-level effects we cannot 
find out whether the charged Higgs emerging from that decay is supersymmetric or not. 
Quantum effects, however, can. 



3. Renormalization of the t6i/+-vertex 

Proceeding closely in parallel with our supersymmetric approach to the conventional 
decay t — > b |1^, [1^, we shall address the calculation of the one- loop corrections to the 
partial width of t — b in the MSSM within the context of the on-shell renormalization 
framework Again we may use both the a or the Gp parametrizations. At one- loop 



order, we shall call the former the "a-scheme" and the latter the "Gi^-scheme". In the 
"a-scheme", the structure constant a = aemil'^ = 0) and the masses of the gauge bosons, 
fermions and scalars are the renormalized parameters: {a, M^, Mz, MH,mf, Msusyi ••■) 
-MsusY standing for the collection of renormalized sparticle masses. Similarly, the "Gi?- 
scheme" is characterized by the set of inputs {Gp, Mw, Mz, MH,mf, Msusy, •••)• Beyond 
lowest order, the relation between the two on-shell schemes is no longer given by eq.(|D 
but by 



V2 2M2 



^For a comprehensive exposition, see e.g. Refs.[p9|-[[4l| 
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where Ar*^'^'^^^ is the prediction of the parameter Ar |^ in the MSSMQ. 

Let us sketch the renormahzation procedure affecting the parameters and fields related 
to the 1 6 if^-vertex, whose interaction Lagrangian was given on eq.(|^. In general, the 
renormalized MSSM Lagrangian C —>■ C + 6C is obtained following a similar pattern as in 
the SM, i.e. by attaching multiplicative renormahzation constants to each free parameter 

1 /2 

and field: gi — > (1 + Sgi/gi)gi, $j Z^^ $j. As a matter of fact, field renormahzation 
(and so Green's functions renormahzation) is unessential and can be either omitted or be 
carried out in many different ways without altering physical (^-matrix) amplitudes. In 



our case, in the line of Refs. [p!7| , |18[ we shall use minimal field renormahzation, i.e. one 
renormahzation constant per gauge symmetry multiplet . In this way the counterterm 
Lagrangian, 6C, as well as the various Green's functions are automatically gauge- invariant. 
Specifically, for the quark fields under consideration, we have 




tR 



L) 
tR. 



(19) 



Here Zi = 1+SZi are the doublet (Zl) and singlet (Z^J^) field renormahzation constants for 
the top and bottom quarks. Although in the minimal field renormahzation scheme there is 
only one fundamental constant, Zl, per matter doublet, it is useful to work with Z'l = Zl 
and Zj^, where the latter differs from the former by a finite renormahzation effect |^0|- To 



fix all these constants one starts from the usual on-shell mass renormahzation condition 
for fermions, /, together with the "residue = 1" condition on the renormalized propagator. 
These are completely standard procedures, and in this way one obtains 



Sm4 



nit 



and 



(20) 



S{,^(mJ)+mJ[S{'(mJ) + 4' 



mi)] 



(21) 



where we have decomposed the fermion self-energy according to 

S^(p) = S{(p2) ^P^ + si(/) ^Pn + nif Eiip') , (22) 

and used the notation S'(p) = dT,{p)/dp'^. 

''A dedicated study of Ar'^'^^^^' has been presented in Ref . 

^We follow the notation of Ref. ||l7|] , which is close enough to that of Ref.|^^, but differs from it 
in several respects, in particular in the sign conventions for the self-energy functions. Moreover, we 
understand that in all formulas defining counterterms we are taking the real part of the corresponding 
functions. 
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One also assigns doublet renormalization constants to the two Higgs doublets (|^) of 
the MSSM: 




Hi 
Hi 



yl/2 

'H2 



HI 



(23) 



The renormalization of the gauge sector is related to that of the Higgs sector. In particular, 
we point out the presence in our decay process t H^ b of the (one-loop induced) mixing 
term H^ — for the bare fields (Cf. Appendix B), which must be renormalized away 
for the physical fields H^ and W^. In order to generate the corresponding Lagrangian 
counterterm we write 



{Zff'w^±i^-^d..H^ 



w 



Therefore, from 



we obtain 



5C 



HW 



^Wht - 

-iSZnw 



(24) 



(25) 



9 



V2Mw 
9 



d.H-b^'^PLt + h.c. 



y2M, 



H 



w 



mtb Put - nibb Pit 



+ h.c. 



(26) 



and in this way it adopts the form of the original vertex (^. In the above expression 

= 1 + 5ZY is the usual SU (2) ^ gauge triplet renormalization constant given by the 
formula 

S^(A;2| 

Sw 



5ZY 



2 
W 



fc2=0 



Ml 



■'W 



'5 Ml 



6M^ ' 
Mir , 



(27) 



and 



5Ml, 



w) 



SMl 



-T.zik'^ = Ml 



(2^ 



are the gauge boson mass counterterms enforced by the usual on-shell mass renormaliza- 
tion conditions. Furthermore, 5Zhw on eqs.(P^-(PBD is a dimensionless constant associ- 
ated to the wave-function renormalization mixing among the bare if^ and fields. Its 
relation with the doublet renormalization constants, ZHi = 1 + 5^/^^., is the following (see 
Appendix B): 

(5 tan/? 



5Zhw = sin P cos P 



2 (^Zh, 



SZh^ + 



tan (3 



(29) 



where 5 tan /3 is a counterterm associated to the renormalization of tan (see below) . 

In practice, the most straightforward way to compute SZhw is from the unrenor- 
mahzed mixed self-energy T^nwik'^) in the unitary gauge, where it takes on the simplest 
form: 

SZ,^ = ^I^eIMM . (30) 
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However, since we shall perform the rest of the calculation in the Feynman gauge , it 
is worth considering the computation of SZh]^ in that gauge (see Appendix B), where 
the discussion is slightly more complicated due to the presence of Goldstone bosons (G^) 
leading to additional (if^ — G^) mixing terms among the bare fields. The corresponding 
expression for 6Zhw is, however, formally identical in both gauges. 
For the SU{2)i gauge coupling constant, we have 

9^{l+^j)9 = {ZY){Zr)-"'9. (31) 

where refers to the renormalization constant associated to the triple vector boson 
vertex. Therefore, from charge renormalization, 



8a _ E^ik"^ 



a k"^ 



-2^21%^, (32) 



and the bare relation a = / At: ^ a + 5a = {g^ + 5g^) (s^ + 5s^)/47r, one gets for 
the counterterm to g: 

bg" 5a cl, (5Ml 5Ml\ 



g^ a si, \ Ml Ml 



(33) 



w 



and clS cl by-product 



5ZY = -^-^ + -5ZY . (34) 



Let us now outline the renormalization of the Higgs sector of the MSSMp5|, De- 
pending on the particular problem at hand, the renormalization procedure may adopt the 
CP-odd state as the basic field on which to set the mass and wave-function renor- 
malization conditions. In the present work, however, since the external Higgs particle is 
charged, we rather take as the basic field. Its mass and field renormalization constants 
are defined by 

M^± -> M^± + 5Mli , ^ Z]ilH^ . (35) 

The charged Higgs field renormalization constant, Zh± = 1 + 5Zh±, is of course related 
to the fundamental doublet renormalization constants introduced on eq.(^), as follows 
(Cf. Appendix B): 

5Zh± = sin^ P 5Zhi + cos^ P 5Zh2 ■ (36) 
The structure of the renormalized self-energy is 

S^±(A;2) = Sh±(A;2) + 5M|± - (P - M|.±) 5Zh± , (37) 

where T,H±{k'^) is the corresponding unrenormalized self-energy. 

In order to determine the counterterms, we impose the following renormalization con- 
ditions: 
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i) On-shell mass renormalization condition: 

Eh±(M|±) = 0, (38) 

ii) "Residue = 1" condition for the renormalized propagator at the pole mass: 



From these conditions one derives 



^S'^4Mi±)=0. (39) 



k2=M^ . 



5Mi± = -S^±(Mi±), 

6Zi,± = +S'^±(M^±). (40) 

Although not needed in our calculation, it is clear that with these settings the neutral 
Higgs fields will undergo an additional finite wave function renormalization. 



Consider next the renormalization of the Higgs potential in the MSSM, eq.(P |25|1 . 
After expanding the neutral components and if^ around their VEV's Vi and V2, 
the one-point functions of the resulting CP-even fields are required to vanish, i.e. the 
tadpole counterterms are constrained to exactly cancel the tadpole diagrams, so that 
the renormalized tadpoles are zero and the quantities vi^2 remain as the VEV's of the 
renormalized Higgs potential. Notwithstanding, at this stage a prescription to renormalize 
tan/3 = V2/V1, 

tan /? — i> tan /3 -|- (5 tan /? , (41) 

is still called for. There are many possible strategies. The ambiguity is related to the 
fact that this parameter is just a Lagrangian parameter and as such it is not a physical 
observable. Its value beyond the tree- level is renormalization scheme dependent. (The 
situation is similar to the definition of the weak mixing angle 6w, or equivalently of 
sin^ 9w-) However, even within a given scheme, e.g. the on-shell renormalization scheme, 
there are some ambiguities that must be fixed. For example, we may wish to define tan /3 
in a process-independent ("universal") way as the ratio V2/V1 between the true VEV's 
after renormalization of the Higgs potential 43|. In this case a consistent choice (i.e. 



a choice capable of renormalizing away the tadpole contributions) is to simultaneously 
shift the VEV's and the mass parameters of the Higgs potential, eq.(|), 

ml2 ^ 4,4,(^?2 + 5^i2), (42) 
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(z = 1,2) in such a way that 5vi/vi = 6v2/v2- This choice generates the following 
counterterm for tan/3 in that scheme (Cf. Appendix B): 

Nevertheless, this procedure looks very formal and one may eventually like to fix the 
on-shell renormalization condition on tan in a more physical way, i.e. by relating it 
to some concrete physical observable, so that it is the measured value of this observable 
that is taken as an input rather than the VEV's of the Higgs potential. Following this 
practical attitude, we choose as a physical observable the decay width of the charged 
Higgs boson into r-lepton and associated neutrino: r+z/^. As it has been argued in 

Section 2, this should be a good choice, because: i) When t if + b is allowed, the decay 
T^Vr is the dominant decay of already for tan f3 ^ 2; ii) From the experimental 
point of view there is a well-defined method to separate the final state r's originating 
from iJ"^-decay from those coming out of the conventional decay — t^i't, so that 
T^Vj. should be physically accessible; and iii) At high tan/5, the charged Higgs 
decay of the top quark can have a sizeable branching ratio. 

The interaction Lagrangian describing the decay r^Vr is directly proportional 

to tan /3, 

q itLt tan /3 , , , 

and the relevant decay width is proportional to tan^ j3. Whether in the a-scheme or in 
the G'iT'-scheme, it reads: 

^<^^ - ^^-) - '-^^ - '^'^ '-^"'^ - ^^"''"> • <^^> 

where we have used the relation (|l^). By measuring this decay width one obtains a physi- 
cal definition of tan (3 which can be used beyond the tree-level. A combined measurement 
of Mh± and tan/? from charged Higgs decaying into r-lepton in a hadron collider has 
been described in the literature |]4^ , ^ by comparing the size of the various signals for 
charged Higgs boson production, such as the multijet channels accompanied by a r-jet or 
a large missing pt, and the two- r-jet channel. At the upgraded Tevatron, the conventional 
mechanisms gg{qq) ti followed by t ^ b have been studied and compared with the 
usual t b, and the result is that for Mjj± ~ 100 GeV the charged Higgs production 

is at least as large as the production, apart from a gap around tan /3 ~ 6 |Q . 

Insofar as the determination of the counterterm 6 tan (3 in our scheme, it can be fixed 
unambiguously from our Lagrangian definition of tan j3 on eq. (|4^) and the renormalization 
procedure described above. It is straightforward to find: 

6 tan (3 6v 1 ^ . , 

- — ^ = -^Zh± + cot 13 5Zhw + A, . 46 

tan [j V 2 
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Notice the appearance of the vacuum counterterm 

which is associated to v'^ = vf + f|, and whose structure is fixed from eq.(^. The last 
term on eq. (|46|) , 

A. = -^-l<5Zr-^5Z]j-F., (48) 

is the (finite) process-dependent part of the counterterm. Here 6mr/mr, 6Z'^£' and are 
obtained from eqs.(pOl) and (with = ); they represent the contribution from 
the mass and wave-function renormahzation of the (z/^, r)-doublet, including the finite 
renormalization of the neutrino leg. Finally, on eq.(^8D is the form factor describing 
the vertex corrections to the amplitude of r^Vr- 

On comparing eqs.(^3D and (|46|) we see that the first definition of tan/3 appears as 
though it is free from pro cess- dependent contributions. In practice, however, process- 
dependent terms are inevitable, irrespective of the definition of tan/3. In fact, the defini- 



tion of tan/? where 6vi/vi = 6V2/V2 |^ will also develop pro cess- dependent contributions, 
as can be seen by trying to relate the "universal" value of tan (3 in that scheme with a 
physical quantity directly read off some physical observable. For instance, if M^o is heavy 
enough, one may define tan (3 as follows: 



l^" mj:. 



r(AO ^ tt) m\ \ M 



V ^2 , 



(49) 



where we have neglected ^ M^o, and bV stands for the vertex corrections to the 
decay processes ^ hh and ^ tt. Since the sum of the mass and wave-function 
renormalization terms along with the vertex corrections is UV-finite, one can consistently 
choose bvxjvx = 6V2/V2 leading to eq. (|43|) . Hence, deriving tan/3 from eq.(|49|) unavoidably 
incorporates also some process-dependent contributions. 

Any definition of tan/3 is in principle as good as any other; and in spite of the fact 
that the corrections themselves may show some dependence on the choice of the particular 
definition, the physical observables should not depend at all on that choice. However, it 
can be a practical matter what definition to use in a given situation. For example, our 
definition of tan /3 given on eq.(|^) should be most adequate for Mh± < mt — rrib and large 
tan/3, since then r+ Ur is the dominant decay of H^, whereas the definition based 

on eq. (|49|) requires also a large value of tan j3 (to avoid an impractical suppression of the 
bb mode); moreover, in order to be operative, it also requires a much heavier charged 
Higgs boson, since M^± ~ M^o > 2mt when the decay A ^ tt is kinematically open in 
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the MSSM. (Use of light quark final states would, of course, be extremely difficult from 
the practical point of view.) 

Within our context, we use eq.(^6D for 6ta.n/3/ tan (3 in order to compute the one-loop 
corrections to our decay t —>■ b. Putting all the pieces together, the counterterm 
Lagrangian for the vertex tbH'^ follows right away from the bare Lagrangian (|^ after 
re-expressing everything in terms of renormalized parameters and fields in the on-shell 
scheme. It takes on the form : 

SCmt = J, H- b [SCr nit cot {3 Pr + 5Cl m, tan /5 P^] t + h.c. , (50) 

V 2 My/ 

with 

= '-^-'-^+'-6Z,.+ '-6Zi+'-6ZU'-^-SZ,^cotp, (51) 
777-^ V z z z tan p 

and where we have set Vtb = 1 [Vtb = 0.999 within ±0.1%, from unitarity of the CKM- 
matrix under the assumption of three generations). 

4. One-loop Corrected r{t H+ b) in the MSSM 

As stated in Section 2, the study of the decay t b is worthwhile in the small 

(tan/3 < 2), and most conspicuously in the high (tan/5 > 30) tan/? region, where the 
branching ratio can be comparable to the one of the standard decay t b. These 

are, therefore, the regions on which we will focus our search for potentially significant 
(strong and electroweak like) SUSY quantum corrections to t — > b. As for the strong 



effects, they can be rather large and have been evaluated in Ref. [|T5|; here we shall not 
dwell any longer on their detailed structure apart from including them in our numerical 
analysis and adding some useful remarks in Section 5. 

On the electroweak side, one may also expect sizeable quantum corrections from en- 
hanced Yukawa couplings of the type (0). In the relevant tan/? regions mentioned above, 
the latters yield the leading electroweak contributions and in these conditions we will 
neglect the pure gauge corrections from transversal gauge bosons in the Feynman gauge. 
Moreover, as already stressed in Section 2, the branching ratio of the charged Higgs mode 
in the intermediate tan (3 region is too small to speak of, so that the detailed structure of 
the radiative corrections in this range is irrelevant. 

In the following we will describe the relevant electroweak one-loop supersymmetric 
diagrams entering the amplitude of t ^ b in the MSSM. At the tree-level, the only 
Feynman diagram is the one in Fig.l. At the one- loop, we have the diagrams exhibited in 
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Figs. 2-6. The computation of the one-loop diagrams requires to use the full structure of 
the MSSM Lagrangian. The explicit form of the most relevant pieces of this Lagrangian, 
together with the necessary SUSY notation, is provided in Appendix A. 

Specifically, Fig. 2 shows the electroweak SUSY vertices involving squarks, charginos 
and neutralinos. In all these diagrams a sum over all indices is taken for granted. The 
super symmetric Higgs particles of the MSSM and Goldstone bosons (in the Feynman 
gauge) contribute a host of one-loop vertices as well (see Fig.3). As for the various self- 
energies, they will be treated as counterterms to the vertices. Their structure is dictated 
by the Lagrangian (|50|) . Thus, Fig.4 displays the counterterms Chi — Ct4 generated from 
the external bottom and top quark lines; they include contributions from supersymmetric 
particles, Higgs bosons and Goldstone bosons. Similarly, Fig. 5 contains the counterterms 
Chi — Cha associated to the self-energy of the external charged Higgs boson. A variant 
of the latter contribution is the mixed — self-energy counterterms Cmi — Cms 
shown in Fig. 6. 

Although we have displayed only the process dependent diagrams, the full analy- 
sis should also include the SUSY and Higgs/Goldstone boson contributions to the vari- 
ous universal vacuum polarization effects comprised in our counterterms. However, the 
calculation of all these pieces has already been discussed in detail long ago in the lit- 
erature [|^, ^ and thus the lengthy formulae accounting for these results will not be 
explicitly quoted here. Their contribution is not tan /^-enhanced, but since we wish to 
compute the full supersymmetric contribution in the relevant regions of the MSSM pa- 
rameter space, those effects will be included in our numerical code. Finally, the smaller 
-though numerically overwhelming - subset of strong supersymmetric one-loop graphs 
are displayed in Fig. 2 of Ref.|]TB[. We will use the formulae from the latter reference in 
the present analysis to produce the total (electroweak -l-strong) SUSY correction to our 
process. 

Next let us report on the contributions from the various vertex diagrams and countert- 
erms in the on-shell renormalization scheme. The generic structure of any renormalized 
vertex function. A, in Figs. 2-3 is composed of two form factors Fl, plus the countert- 
erms. Therefore, on making use of the formulae of Section 3, one immediately finds: 

A = [rrit cot/5 (1 + Kr) Pr + tan/5 (1 + A^) , (52) 

V 2 Mw 



where 



= FR+^-^+'-SZi+'-SZ^^-A^ 
— —5- + SZh+ + (tan f3 — cot P) SZ^w 

V 
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Al = Fi + ^ + ^5Zi + ^54 + A.. (53) 

In the following the analytical contributions to the vertex form factors and counterterms 
will be specified diagram by diagram. 

4.1 SUSY vertex diagrams 

In this section we will make intensive use of the definitions and formulae of Appendix A. 
We refer the reader there for questions about notation and conventions. Following the 
labelling of Feynman graphs in Fig. 2 we write down the terms coming from virtual SUSY 
particles. 



Diagram (V51): Making use of the coupling matrices of eqs. (|A. 18|) and ( |A.23| ) we 
introduce the shorthand^ 

A^^A^ll and ^ 4L , (54) 

and define the combinations (omitting indices also for Qai^Qai) 

= cos pAlQ^A^^^ , ^(1) = cos pA*_Q^A^^^ , 
B^'^ = cosPAlQ^Af , F(^) = cos pA*^Q^Af , 

c(^) = sin/?A;g^AL°^ , = sin/5A*_g^AL°) , 

D(^) = sin (3 AlQ^Af , H^^'^ = sin (3 A*_Q^Af . (55) 

The contribution from diagram (Ksi) to the form factors Fl and Fr is then 

Fl = Ml[h^'^Co+ 

+ mt [mtH^^^ + M°G(i) + m^A^^^ + MiE^^^) (Cn - C12) 
+ fm^m^A^^) + mtMiE^^^ + M^rubB^^^ + MiM^F^^A Co 



= Mr 



+ rrik {mtH^^^ + M^G^^^ + mfeA^^^ + MiE^^^) C12 

+ mt (mtA(i) + M05(i) + mfeiy(i) + M^D^i)) (Cn - C12) 

[mtrribH^^^ + mtMiO'^^^ + Mlm^G^^^ + MiM^C^^A Cq] , (56) 



where the overall coefficients Ml and M/j are the following: 



ig'^Mw ig'^M^ 



Ml = -^^^^-^ Mr = -^^^. (57) 
■rrib tan p rrit cot p 

^Lower indices are summed over, whereas upper indices (some of them within parenthesis) are just 
for notational convenience. 
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The notation for the various 3-point functions is summarized in Appendix C. On 
eq. (|56|) they must be evaluated with arguments: 



(58) 



Diagram (V52): For this diagram -which in contrast to the others is finite- we also 
use the matrices on eqs. ( |A.18| ) and ( |A.21| ), and introduce the shorthands 



A 



Ai'L and A 



(t) 



A 



(59) 



to define the products of coupling matrices 



= GtaA^t^*A^l^ , = GkaA^^^*A^l^ , 

5(2) = G,aA?*A^l^ , = . 

The contribution to the form factors and from this diagram is 

m,5(2)Ci2 + miC(2) (Cn - C12) - M^^D^'^Co 
m,C^^^Ci2 + (Cn - C12) - M^A^^^Co 



(60) 





Ml r 


2Mw ^ 


Fr = 


Mr r 


2Mw. ^ 



(61) 



the coefficients iW^L, Mr being those of eq. (0) and the scalar 3-point functions now 
evaluated with arguments 

(62) 



(63) 



= [p, p,M^, mi^ , mi^ 

Diagram (V53): For this diagram we will need 

A± = A^^\^ and A^^.^ = A±^aa ^ 
and again omitting indices we shall use 

= cosM°)*gM_, 

=cos/5Af*Q^A+, 
C^'^ = sin PA^^^*Q^ A., 

From these definitions the contribution of diagram (¥33) to the form factors can be 
obtained by performing the following changes in that of diagram (V51), eq. (^6]): 

— Everywhere on eqs. (^) and (|58D replace Mj <-> M° and m^^ ^ m^^. 

— Replace on eq. (|56D couplings from (|55D with those of (|6^) . 



^(3) = cosM°^*gM_, 

=cosM°^*Q^A+, 
G^^') =sinM°)*Q^A_, 



(64) 
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— Include a global minus sign. 



4.2 Higgs vertex diagrams 

Now we consider contributions arising from the exchange of virtual Higgs particles and 
Goldstone bosons in the Feynman gauge, as shown in Fig.3. We follow the vertex formula 
for the form factors by the value of the overall coefficient N and by the arguments of the 
corresponding 3-point functions. 

• Diagram (V^i): 

Fl = 7VK(Ci2-Co)+m?cot2/?(Cn-Ci2)], 

Fr = Nml[Cu-Co + tan^P{Cn-Cu)], 

ig^ A _ {M|o^M|o}\ {cos g, sin g} 
2 \ 2M^ ) cos/5 



Diagram (Vf/2): 

Fl = 7Vcot/3K(Cii-Ci2)+m^(Co-Ci2)], 
Fr = 7Vm^tan/3(2Ci2-Cii-Co), 

= a (p, p', mfe, M^., {Mjfo , Mfto }) . 

Diagram (Vffs): 

Fl = iVm,2[cot=^/3Ci2 + Cn-Ci2-Co], 

F« = 7V[m2tan2/3Ci2 + m2(Cn-C/i2-Co)], 

igf^ {sina, cosa} \ ■ //p m /i {^^0'^fto}\ 



2 



= - — ^{cos(/3-Q;),sin(/3-Q;)} 1 

2 sin/3 ^ ^' ^-"V 2M^^ 

= a(p,p',mt,{Mj^o,M^o},M^±) . 



Diagram (Vh4): 

Fl = Arm?(2Ci2-Cii + Co)cot/3, 

= iV[-m2Ci2 + m,2(Cn-Ci2-Co)]tan/5, 
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Diagram (V^s): 



Fl = iVK(Ci2 + Co)+m,^(Cn-Ci2)] 
Fr = iVm^an2/3(Cn + Co) , 



N 



AO 



4 M2,^ 



= C* (p, p', TUb, Mw, Mao) . 



Diagram (V^g) 



Fl = Nm^,coef3{Cu + Co), 

Fr = N[mlCi2 + ml{Cii-Ci2 + C^)] 



N 



ig" (Ml^ M\; 
4 \Ml, M^^ 



= {p, p',mt, Mao , Mw) ■ 



Diagram {Vm)'- 



Fl = N[{2mlCu + Co + 2im^,-ml)iC^,-Cu))cot^P + 2ml{Cu + 2Co)]ml 

Fr = N[{2mlCn + Co + 2{ml-ml){Cn-C,2))t^n^(3 + 2ml{Cu + 2Co)]ml, 

, iq"^ sin a cos a 

N = ± — 

4M^ sin /3 cos /3 ' 

= C^{j),p\{MHO,Mho},mt,mb) . 



Diagram {Vm) 



Fl = Nm] coi^ (3 Cq, 



Fr = Nmltan'^pCo, 
N = t''^ 



= C^{p,p',{MAo,Mz},mt,mb) . 
In the equations above, it is understood that the CP-even mixing angle, a, is renor- 



mahzed into aes by the one-loop Higgs mass relations p5 . 

As for the SUSY and Higgs contributions to the counterterms, they are much simpler 
since they just involve 2-point functions. Thus we shall present the full electroweak results 
by adding up the various sparticle and Higgs effects. In the following formulae, we append 
labels referring to the specific diagrams on Figs. 4-6. 
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4.3 Counterterms 

• Counterterms 6m f ,6Z[ ,SZ{^: For a given down-like fermion b, and corresponding 
isospin partner t, the fermionic self-energies receive contributions 



A 



it) 



iCbl) + iCb2) 



A 



(b) 
itao 



2 



+ ^M°Re Bo {p, M° 



(65) 



from SUSY particles, and 



{ 



(C63) + {C64) 
2 



{cot'/5, tan'/3}5i(p, m^, Mj^±) + m*, M^) 



+ 



cos^a Bi{p, nib, M^o) + sin^a Bi{p, nib, M^o) 



2 cos2/? 

+ sin2/35i(p,mb,MAo) + cos^/?5i(p,m6,Mz)]} 



(C63) + (Cm) 



+ 



{m^ [Bo{p, mt, Mh±) - Bo{p, mt, Mw)] 

cos^a Bo{p, mb, M^o) + sin^a Bq{p, mb, M; 
2cos''p L 

sin^/3 5o(p, mfc, M^o) - cos^/? 5o(p, m^, M^) | 



mi 



(66) 



from Higgs and Goldstone bosons in the Feynman gauge. To obtain the correspond- 
ing expressions for an up-like fermion, t, just perform the label substitutions b ^ t 
on eqs. (|H3|)-(|0BD; and on eq. (|5^) replace sin a cos a and sin/5 cos/5 (which 
also implies replacing tan/5 cot /?). 

Introducing the above expressions into eqs. ([20|)-(^TD one immediately obtains the 
SUSY contribution to the counterterms 6m f , SZ[ pi^. 

Counterterm 6Zh±: 



5Zfj± = 5Zj{± 



{CHl) + {CH2) + {CH'i) + {CH4) 
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{ml tan^/^ + m] cot^p){Bi + Ml±B[ + mlB'^) 
+ 2mlmfBo {Mh± , m^, m*) 
NcJ2\GbafB',{MH±,m-,^ 

ab 



gSf Sin'/?) (5i + Ml^B[ + MfB',) 



{Mh±,M'^,M,[ 



+ 2MiM^Re (Q^iQ^*) sin/3cos/3E^ 
Notice that diagram (Chs) gives a vanishing contribution to 5Zj:/± 
Counterterm 6Zh\y: 

SZhwI (Cmi )+ (Cm2 )+ (Cm3 ) 



(67) 



5Z 



HW 



M? 



w 



= -^Jr^\mltanP{Bo + Bi) +mlcotpBi {MH±,mi„mt) 
'^'^^T.GbaR^aR'S* [2B, + Bo] {MH.,m~,^,mO 

W ab 



2M? 



+ M, (sin/3 Q^*C^, + cos/3 Q^*^^) 5i] (M^±, M°, M,; 
A sum is understood over all generations. 



(68) 



Finally, the evaluation of A,- on eq.(P5|) yields similar bulky analytical formulae, which 
follow after computing diagrams akin to those in Figs. 2-6 for the MSSM corrections to 
— s> r+ Ut-. We refrain from quoting them explicitly here. The numerical effect, though, 
will be explicitly given in Section 5. 

We are now ready to furnish the corrected width of t ^ b in the MSSM. It just 
follows after computing the interference between the tree-level amplitude and the one-loop 
amplitude. It is convenient to express the result as a relative correction with respect to 
the tree- level width both in the a-scheme and in the G^-scheme. In the former we obtain 
the relative MSSM correction 

r - r(°) 



rMSSM 



p(0) 
i a 



Lr. ■ -«roo... ^. , [2 i?e(Ai + A^)] , (69) 



^[2i?e(A,)] + ^[2i?e(A«)]+ ^ 



where the corresponding lowest-order width is 

D 



r(o) _ ( ^ 



1/2 m| M^±, 



2 ' 2 

rrif mf 



(70) 
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with 



D = {ml + ml- Mh± ) {ml cot^ {3 + ml tan^ j3) + Amlml , 

Nl = {ml + ml- Mh± ) ml tan^ p , 

Nr = {ml + ml- Mh± ) ml cot^ p , 

Nlr = 2mlml. (71) 

From these equations it is obvious that at low tan (3 the relevant quantum effects basically 
come from the contributions to the form factor Ar whereas at high tan (3 they come from 

Using eq.([18D we find that the relative MSSM correction in the G^-parametrization 
reads 

r - r(°) 

;:MSSM _ ^ ^ Gp _ rA/S5M a MSSM (^r,\ 
'^Gf - p(0) - , [(A) 

^ Gf 

where the tree- level width in the G^-scheme, T^^, is given by eq.(P) and is related to 
eq.(|70D through 

rw = rg^ (1 - Ar*^^^*0 • (73) 
5. Numerical Analysis and Discussion 

Quantum effects should be able to discriminate whether the charged Higgs emerging 
from the decay t b is supersymmetric or not, for the MSSM provides a well defined 

prediction of the typical size of these effects using the present bounds on sparticle masses. 
Some work on radiative corrections to the decay width of t ^ H'^ b has already appeared 
in the literature. In particular, the conventional QCD corrections have been evaluated [T^ 



and found to significantly reduce the partial width. The SUSY-QCD corrections are also 
substantial and have been analyzed, only in part in Refs. P8|, and in more detail 
in Ref. ||15||. Nevertheless, the electroweak corrections produced by the roster of genuine 



(R-odd) sparticles have not been considered at all yet. As for the virtual effects medi- 
ated by the Higgs bosons, a first treatment is given in Refs.[5C] and [^]. However, these 
references disagree in several parts of the calculation, and moreover they are both incom- 
plete calculations on their own, for they fully ignore the Higgs effects associated to the 
bottom quark Yukawa coupling, which could in principle be significant in the large tan /3 
region. On the other hand, even though the latter kind of Higgs effects have been dis- 
cussed in the literature in other renormalization schemes based on alternative definitions 
of tan P a detailed analysis including the genuine SUSY effects themselves has 

never been attempted. Thus, if only for completeness, we are providing here not only a 
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dedicated treatment of the fZ-odd contributions mediated by the sparticles of the MSSM, 
but also the fully-fledged pay-off of the supersymmetric Higgs effects. 

Before presenting the results of the complete numerical analysis, it should be clear 
that the bulk of the high tan (3 corrections to the decay rate of t ^ b in the MSSM is 
expected to come from SUSY-QCD. This could already be foreseen from what is known 
in SUSY GUT models |]lT|; in fact, in this context a non- vanishing sbottom mixing (which 
we also assume in our analysis) may lead to important SUSY-QCD quantum effects on 
the bottom mass, mb = mf'^'^ + Amfe, where Arrib is proportional to — /itan/3 

at sufficiently high tan f3. These are finite threshold effects that one has to include when 
matching the SM and MSSM renormalization group equations (RGE) at the effective 
supersymmetric threshold scale, Tsusyi above which the RGE evolve according to the 



MSSM /3-functions in the MS scheme [§4|, In our case, since the bottom mass is an 



input parameter for the on-shell scheme, these effects obviously have a different physical 
meaning, but are formally the same; they are just fed into the mass counterterm 6mb/mb 
on eq.(|53|) and contribute to it with opposite sign {5mb/mb = — Am^ -|- ...) []. 

Explicitly, when viewed in terms of diagrams of the electroweak-eigenstate basis, the 
relevant finite corrections from the bottom mass counterterm are generated by mixed 
LR-sbottoms and gluino loops (Cf. Fig. 7a): 

3^ ^9 ^LR n^br ' ' ^g) 

■rrig utanP I{m-^^,mi^,mg) , (74) 




where the last result holds for sufficiently large tan (3 and for fi not too small as compared 
to Ab. We have introduced the positive-definite function (Cf. Appendix C) 

2 2 1 Tl? I 2 2 1 '"o I 2 2 1 "l-o 

mf mi m—j + mi mi m — | H- mf mi m — f 

J(mi, ma, mg) = 16 tt « Co(0, 0, mi, ma, mg) = — 2\ ^ 2 ow 2 2^ ■ 

(m| — ma) (ma — m^) (mf — m^) 

(75) 

In addition, we could also foresee potentially large (finite) SUSY electroweak effects from 

6mb/mb. They are induced by tan /^-enhanced Yukawa couplings of the type @. Of 

course, these effects have already been fully included in the calculation presented in Section 

3 that we have performed in the mass-eigenstate basis, but it is illustrative of the origin 

■^In the alternative framework of Ref.||l|l, the SUSY-QCD corrections have been computed assuming no 
mixing in the sbottom mass matrix. Nonetheless, the typical size of the SUSY-QCD corrections does not 
change as compared to the present approach (in which we do assume a non-diagonal sbottom matrix) the 
reason being that in the absence of sbottom mixing, i.e. Af]^^ = 0, the contribution Smb/mt, oc — /xtan/3 
at large tan/3 is no longer possible but, in contrast, the vertex correction does precisely inherits this 
dependence and compensates for it (Cf. Appendix A). The drawback of an scenario based on M]^^ = 0, 
however, is that when it is combined with a large value of tan /3 it may lead to a value of Ah which 
overshoots the natural range expected for this parameter. 
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of the leading contributions to pick them up again directly from the diagrams in the 
electroweak-eigenstate basis. In this case, from loops involving mixed LR-stops and mixed 
charged higgsinos (Cf. Fig. 7b), one finds: 

(^) = — Mlj^Iimi^^mi^, n) 

V^''/ SUSY- Yukawa l^VT^ m^, 

* fi tan p At I {mi^,mi^,fj.) , (76) 



167r 

where again the last expression holds for large enough tan f3. 

Notice that, at variance with eg . (|7^ , the Yukawa coupling correction (|76D dies away 
with increasing fi. Setting ht ^ 1 at high tan/3, and assuming that there is no large 



hierarchy between the sparticle masses, the ratio between ([7^ ) and (^) is given, in good 
approximation, by Arrig/At times a slowly varying function of the masses of order 1, where 
the (approximate) proportionality to the gluino mass reflects the very slow decoupling rate 
of the latter [ID. 

In view of the present bounds on the gluino mass, and since At (as well as Ah) can- 
not increase arbitrarily (Cf. Appendix A), we expect that the SUSY-QCD effects can 
be dominant and even overwhelming for sufficiently heavy gluinos. Unfortunately, in 
contradistinction to the SUSY-QCD case, there are also plenty of additional vertex con- 
tributions both from the Higgs sector and from the stop-sbottom/gaugino-higgsino sector 
where those Yukawa couplings enter once again the game. So if one wishes to trace the 
origin of the leading contributions in the electroweak-eigenstate basis, a similar though 
somewhat more involved exercise has to be carried out also for vertex functions. Of course, 
all of these effects are automatically included in our calculation of Section 3 within the 
framework of the mass-eigenstate basisQ. 

We may now pass on to the numerical analysis of the over-all quantum effects. After 
explicit computation of the various loop diagrams, the results are conveniently cast in 
terms of the relative correction with respect to the tree-level width: 



V - rg^ _ r(t ^ H+ h) - r(o) (t ^ h+ h) 



(77) 



In what follows we understand that 5 defined by eq.(^) is 5a - Cf. eq.(^) - i.e. we shall 

always give our corrections with respect to the tree-level width r° in the a-scheme. The 

corresponding correction with respect to the tree-level width in the Gi^-scheme is simply 

given by eq. ([7^) , where Ar^'^^^'^'^ was object of a particular study and therefore it 

^The mass-eigenstate basis is extremely convenient to carry out the numerical analysis, but it does 
not immediately provide a "physical interpretation" of the results. The electroweak-eigenstate basis, in 
contrast, is a better bookkeeping device to trace the origin of the most relevant effects, but as a drawback 
the intricacies of the full analytical calculation can be (in general) abhorrent. 
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can be easily incorporated, if necessary. Notice, however, that Ar*-'^"^'^^^ is already tightly 
bound by the experimental data on Mz = 91. 1863 ± 0.0020 GeV^ at LEP [^] and the ratio 
Mw/Mz in pp, which lead to Mw = 80.356 ± 0.125 GeV. Therefore, even without doing 
the exact theoretical calculation of Ar within the MSSM, we already know from 

Ar = 1 - ^ I (78) 

that Ar'^^^^ must lie in the experimental interval Ar"'''^ ~ 0.040 ± 0.018. 

Now, since the corrections computed in Section 3 can typically be about one order of 
magnitude larger than Ar^^^^, the bulk of the quantum effects on t ^ b is already 
comprised in the relative correction ( |77[ ) in the a-scheme0. Furthermore, in the conditions 
under study, only a small fraction of Ar^^^^^ is supersymmetric [^], and we should not be 
dependent on isolating this universal, relatively small, part of the total SUSY correction 
to 5. To put in a nutshell: if there is to be any hope to measure supersymmetric quantum 
effects on the charged Higgs decay of the top quark, they should better come from the 
potentially large, non-oblique, corrections computed in Section 3. The SUSY effects 
contained in Ar^^^^^'^ , instead, will be measured in a much more efficient way from a 
high precision ((5M^p = ±40 MeV) determination of Mw at LEP 200. 

Another useful quantity is the branching ratio 

^ BRit ^H+b)= ^/""^^ , (79) 

i- W + i- H + i- SUSY 

where Tw = r(i — *■ h) and Tsusy stands for decays of the top quark into SUSY 
particles. In particular, the potentially important SUSY-QCD mode t ^ ti ^ is kinemat- 
ically forbidden in most part of our analysis where we usually assume = (9(300) GeV . 
There may also be the competing electroweak SUSY decays t — > X« aiid t ^ bixf for 
some a = 1, ...,4 and some i = 1,2. The latter, however, is also phase space obstructed 
in most of our explored parameter space, since we typically assume m^^ = 150 GeV. The 
decay t —>■ ti Xa? instead, is almost always open, but it is not tan /9-enhanced in our 
favourite segment (jl^. However, when studying the branching ratio (^) as a function 
of the squark and gluino masses, we do include the effects from all these supersymmetric 
channels whenever they are kinematically open. Thus in general Tsusy on eq. (|79|) is given 

by 

Tsusy = T{t ^ h ^) + E ^ xl) + E r(t ^ bi xt) ■ (80) 



^For the standard decay t b, the situation is quite different since the SM electroweak correc- 

tions 1 57 and the maximal SUSY electroweak corrections [0 in the a-scheme are much smaller than for 



the decay t iJ+ b , namely they are of the order of Ar. Therefore, for the standard decay t — > W'^ b 
there is a significant cancellation between the corrections in the a-scheme and Ar in most of the tan /3 
range resulting in a substantially diminished correction in the Gj^-scheme. 
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The various terms contributing to this equation are computed at the tree-level. Recently, 
the SUSY-QCD corrections to some of these supersymmetric modes have been evaluated 



and in some cases may be important ||5^ . Similarly, we treat the computation of the 
partial width of the standard mode t —>■ b at the tree-level. This is justified since, as 
shown in Refs. O, ITH|, O , this decay cannot in general develop large supersymmetric 



radiative corrections, or at least as large as to be comparable to those affecting the charged 
Higgs mode (for the same value of the input parameters). The reason for it stems from the 
very different structure of the counterterms for both decays; in particular, the standard 
decay mode of the top quark does not involve the mass renormalization counterterms for 
the external fermion lines, and as a consequence the aforementioned large quantum effects 
associated to the bottom quark self-energy at high tan/3 are not possible. 

Figures 8-20 and 22-23 display a clear-cut resume of our numerical results. We wish 
to point out that they have been thoroughly checked. Scale independence of 6, eq. (ffTf) , 
and cancellation of UV-divergences have been explicitly verified. Most of the analytical 
and numerical calculations have been doubled. In particular, we have constructed two 
independent numerical codes and checked that the two approaches perfectly agree at 
different stages. In all our numerical evaluations we have imposed the following restriction 
on the non-SM contributions to the p-parameter 0: 

(5pncw < 0.003 . (81) 

To start with, we concentrate on the case /i < 0, which we study in Figs. 8-20. (The 
case /i > is studied apart in Figs. 22-23 and will be commented later on.) We observe 
that, for negative /i, the leading SUSY-QCD effects on 6 are positive. This means that 
in these circumstances the potentially large strong supersymmetric effects are in frank 
competition with the conventional QCD corrections, which are also very large and stay 
always negative as will be discussed later on. 

Needless to say, a crucial parameter to be investigated is tan/5. In Fig. 8 we plot 
the tree-level width, Toit H^b), and the total partial width, TMssAiit H^b), 
comprising all the MSSM effects, as a function of tan (3. A typical set of parameters is 
chosen well within canonical expectations (see below); the individual influence of each 
one of them is tested in Figs. 10 to 20. Also shown in Fig. 8 is the (tree-level) partial 
width of the standard top quark decay t b, which is (as noted above) far less 

sensitive to quantum corrections. For convenience, we have included in Fig. 8 a plot of 
rQCD(i!^ b), i.e. the partial width that would be obtained in the presence of only the 

standard QCD corrections. In practice, this is tantamount to saying that rQCD(^ —>■ b) 
is the partial width that would be expected in the absence of SUSY effects, for the 
electroweak non-supersymmetric corrections turn out to be negligible versus the ordinary 
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QCD effects. 

From eq.(^) it is clear that, for large (resp. small) tan /5, the renormalized form factor 
yielding the bulk of the SUSY contribution is (resp. Ar). To appraise the relative 
importance of the various types of MSSM effects on r{t —>■ b), in Figs.9a-9b we provide 
plots for the correction to the partial width, eg. ([77D , and to the branching ratio, eg. (|70D , 
as a function of tan/5, reflecting the various individual contributions. Specifically, we 
show in Fig. 9a: 

• (i) The super symmetric electroweak contribution from genuine (i?-odd) sparticles 
(denoted 5susy-ew), i-e. from sfermions (sguarks and sleptons), charginos and 
neutralinos; 

• (ii) The electroweak contribution from non-supersymmetric (i?-even) particles {6ew)- 
It is composed of two distinct types of effects, namely, those from Higgs and Gold- 
stone bosons (collectively called "Higgs" contribution, and denoted ^Higgs) plus the 
leading SM effects from conventional fermions (^sm): 

Sew = (^Higgs + ^SM ; (82) 

The remaining non-supersymmetric electroweak effects are subleading and are ne- 
glected. 

• (iii) The strong supersymmetric contribution (denoted by 5susy-qcd) from sguarks 
and gluinos; 

• (iv) The strong contribution from conventional guarks and gluons (labelled 5qcd); 
and 

• (v) The total MSSM contribution, (5mssm, namely, the net sum of all the previous 
contributions: 

(^MSSM = ^SUSY-EW + ^EW + (^SUSY-QCD + ^QCD- (83) 

In Fig. 9b we reflect the impact of the MSSM on the branching ratio, as a function of 
tan/3; also shown are the tree-level value of the branching ratio and the latter guantity 
after including the (non-supersymmetric) QCD corrections. A typical common set of 
inputs has been chosen in Figs.9a-9b such that the supersymmetric electroweak corrections 
reinforce the strong supersymmetric effects (SUSY-QCD). For this set of inputs, the total 
MSSM correction to the partial width of t ^ b is positive for tan/5 > 20 (approx.). 
Remarkably enough, this is so in spite of the huge negative effects induced by QCD. In 
fact, we see that the gluon effects are overridden by the gluino effects provided tan/5 is 
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sufficiently large, to be concrete for tan/3 > 30. Beyond this value, the strenght of the 
super symmetric loops becomes rapidly overwhelming; e.g. at the representative value 
tan/3 = mt/nif, = 35 we find 5mssm — +27%; and at tan/3 ~ 50, which is the preferred 



value claimed by 50(10) Yukawa coupling unification models [0], the correction is already 
i^MSSM — +55%. Quite in contrast, at that tan P one would expect, in the absence of SUSY 
effects, a (QCD) correction of about —57%, i.e. virtually of the same size but opposite in 
sign!. 

Coming back to Fig.8, we see that, after including the SUSY effects, the partial width 
of t ^ b equals the partial width of the standard decay t b near the "5*0(10)" 

point tan (3 = 50. (The meeting point is actually a bit earlier in tan (3, after taking into 
account the known |1T7|, |T8|, negative, SUSY corrections to t — > b, but this effect is 
not shown in the figure since it is relatively small.) Now, for the typical set of parameter 
values introduced in Fig.8, the top quark decay width into SUSY particles, eq. (^0]) , is 
rather tiny. Thus it is not surprising that in these conditions the branching ratio of 
the charged Higgs mode can be remarkably high: BR{t — » b) ~ 50%, i.e. basically 
50% — 50% versus the standard decay mode. In contrast, the branching ratio without 
SUSY effects (i.e. essentially the QCD-corrected branching ratio) is much smaller: at 
the characteristic 50(10) value, tan/3 = 50, it barely reaches 20%. Clearly, if the SUSY 
quantum effects are there, they could hardly be missed!. 

As noted before, even though the dominant MSSM effects are, by far, the QCD and 
SUSY-QCD ones, they have opposite signs. Therefore, there is a crossover point of the 
two strongly interacting dynamics, where the conventional QCD loops are fully counter- 
balanced by the SUSY-QCD loops. This leads to a funny situation, namely, that at the 
vicinity of that point the total MSSM correction is given by just the subleading, albeit non- 
negligible, electroweak supersymmetric contribution: 5mssm — ^susy-ew- The crossover 
point occurs at tan/3 > 32 ~ rrit/mb, where 5susy-ew ~ 20. For larger and larger tan/3 
beyond mt/rrib, the total (and positive) MSSM correction grows very fast, as we have said, 
since the SUSY-QCD loops largely overcompensate the standard QCD corrections. As a 
result, the net effect on the partial width appears to be opposite in sign to what might 
naively be "expected" (i.e. the QCD sign). Of course, this is not a general result since 
it depends on the actual values of the MSSM parameters. In the following we wish to 
explore the various parameter dependences and in particular we want to assess whether 
a favourable situation as the one just described is likely to happen in an ample portion of 
the MSSM parameter space. In particular, the value tan/3 = mt/rrib = 35 will be chosen 
in all our plots where that parameter must be fixed. We consider it as representative of 
the low end of the high tan/3 segment, eq.(|T^). Thus tan/3 = nit /nib = 35 behaves as a 
sort of threshold point beyond which the MSSM quantum effects on t b take off so 
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fast that they should have indehble experimental consequences on top quark physics. 

As regards to the non-supersymmetric electroweak corrections, ^ew; it is apparent 
from Fig.9a that they are very small, especially in the high tan/5 segment. Also in the 
very low tan/3 segment, 0.5 < tan/3 < 1, 6ew is relatively small; and this is so not only 
because both ^Higgs and 6sm become never too large in absolute value, but also because 
in that region there is a cancellation between ^Higgs < and 6sm > 0. As it happens, we 
end up with the fact that the complicated Higgs effects result in a very tiny contribution, 
except in the very low tan (3 end, where e.g. they can reach —15% at tan/3 ~ 0.5. In this 
corner of the parameter space, ^Higgs becomes the dominant part of 5mssm, being even 
larger than the QCD effects, which stay at the level of —8%, and also larger than the 
SUSY-QCD and SUSY-EW corrections, which remain below +4% and —1%, respectively. 

We have treated in detail the very low tan (3 segment by including the one-loop renor- 



malization of the Higgs masses [^]. This is necessary in order to avoid that the lightest 
CP-even Higgs mass either vanishes at tan /3 = 1 or becomes lighter than the phenomeno- 
logical bounds near that value. In passing, we have checked that the one-loop shift of the 
masses, as well as of the CP-even mixing angle, a, has little impact on the partial width 
of t — b in the entire range of tan/3, eq.(|T3|). They entail at most an additional 5% 
negative shift of 5mssm in the very low tan (3 region]^. It is precisely in this region where 
the Higgs effects could be expected of some relevance, and thus where the renormalization 
of the CP-even mixing angle could have introduced some noticeable change in the neutral 
Higgs couplings. Quite on the contrary, at high tan /3 the corresponding effect is found 
to be of order one per mil and is thus negligibly small. On the other hand, a simple 
inspection of Figs. 8 and 9b shows that even in the very low tan/3 ballpark, where there 
may be some ten percent effect from the Higgs sector, the rising of the tree-level width is 
so fast that it becomes very hard to isolate these corrections. We conclude that, despite 
the rather large number of diagrams involved, the over-all yield from the Higgs sector of 
the MSSM on t ^ b is rather meagre in the whole tan j3 range (|13|). This fact is some- 
what surprising and was not obvious a priori, due to the presence of enhanced Yukawa 
couplings (0) in the whole plethora of Higgs diagrams. The cancellations involved are 
reminiscent of the scanty SUSY Higgs effects obtained for the standard top quark decay 



t W+b^]. 

We come now to briefly discuss the standard QCD effects up to 0{as), which involve 

one- loop gluon corrections and gluon bremsstrahlung fl^. As it is plain from Fig. 9a, 5qcd 

is negative-definite and very important in the high tan/3 segment. It quickly saturates 

^°To perform that check, we have included both the stop and sbottom contributions to the one-loop 
Higgs mass relations. A set of 7 independent parameters has been used to fully characterize these effects, 
viz. (M//± , /i, tan/3, TO^^ , TOjj^ , ^fc, At). We refrain from writing out the cumbersome formulae [p5[ . 
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for tan/5 ^ 10 at a large value of order —60%. Therefore, the QCD effects need to be 



considered in order to isolate the virtual SUSY signature [|14|. The leading behaviour of 
the standard QCD component in the relative correction ( [77| ) can be easily assessed by 
considering the following asymptotic formula 

o ^ ^^(^2 tan2 + cot^ f3) + 3(4 + tan^ f3 - 2^ cot^ (3)ml In ( ^] 

Oqcd = m2tan2/? + mfcot2/3 ' 

(84) 

which we have obtained by expanding the exact one- loop formula up to 0{ml / mf, Mjj+ / mf 
Here as = as(mj), normalized as q;s(M|) ^ 0.12. The big log factor ln(m^/m^) originates 
from the running b-quark mass evaluated at the top quark scale. The correction is seen 
to be always negative. We point out that while we have used the exact 0{as) formula for 
the numerical evaluation, the approximate expression given above is sufficiently accurate 
to convey the general features to be expected both at low and at high tan /3. In particular, 
for rrih 7^ and tan f] in the relevant high segment (P^), the QCD correction becomes very 
large and saturates at the value 



2a, /87r2-15 , 



OQCD = 1 + In ^ ~ -62% (tan/5 » Jmt/nn, ~ 6) . (85) 

Ti \ 6b ml J * 

(The exact 0{as) formula gives slightly below —60%.) At low values of tan/5, the 
corrections are much smaller, as it follows from the approximate expression (5qcd — 
(— a<(/7r)(87r2 — 15)/18 ~ —12%. We remark that for = the dependence on tan/5 
totally disappears from eq.(Q), so that one would never be able to suspect the large 
contribution (|85|) in the high tan/5 regime. The limit = 0, nevertheless, has been 



considered for the standard QCD corrections in some places of the literature |60|, ^ but, 
as we have seen, it is untenable unless one concentrates on values of tan/5 of order 1, 
in which case the relevance of our decay for SUSY is doomed to oblivion. This situa- 
tion is similar to the one mentioned above concerning the SUSY-QCD corrections in the 
limit = 0, which leads to an scenario totally blind to the outstanding supersymmetric 
quantum effects obtained for mj, 7^ at high tan /5 [|l^ . We stress that in spite of the 
respectable size of the standard QCD effects, they become fast stuck at the saturation 
value (^51), which is independent of tan/5. On the contrary, the SUSY-QCD effects grow 
endlessly with tan /5 and thus rapidly overtake the standard QCD prediction. 

Worth noticing is the evolution of the quantities (|77|) and (^) as a function of the 
gluino mass (Cf. Figs.lOa-lOb). Of course, only the SUSY-QCD component is sensitive to 



mg. Although the SUSY-QCD effects have been object of a particular study in Ref. fTS 
we find it convenient, to ease comparison, to display the corresponding results in the very 
same conditions in which the electroweak supersymmetric corrections are presented. The 
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steep falls in Fig. 10a are associated to the presence of threshold effects occurring at points 



satisfying rrig + m^-^ ~ m^. An analogous situation was observed in Ref . ||T7|, |T8[ for the 
SUSY corrections to the standard top-quark decay. Away from the threshold points, the 
behaviour of Ssusy-qcd is smooth and perfectly consistent with perturbation theory. In 
Fig. 10b, where the branching ratio (|7DD is plotted, the steep falls at the threshold points 
are no longer present since they are compensated for by the simultaneous opening of the 
two-body supersymmetric mode t ^ tig, for rrig < nit — "^fj- 

We emphasize that the relevant gluino mass region for the decay t h is not the 

light gluino region, but the heavy one, the reason being that the important self-energy 
correction mentioned above, eg. (|71D , involves a gluino mass insertion. As a consequence, 
virtually for any set of MSSM parameters, there is a well sustained SUSY correction for 
any gluino mass above a certain value, in our case > 250 — 300 GeV . The correction 
raises with the gluino mass up to a long flat maximum before bending -very gently - 
into the decoupling regime (far beyond ITeV). The fact that the decoupling rate of the 
gluinos appears to be so slow has an obvious phenomenological interest. 

Next we consider in detail the sensitivity of our decay on the higgsino-gaugino param- 
eters (/i, M) characterizing the chargino-neutralino mass matrices (Cf. Appendix A). We 
start with the supersymmetric Higgs mixing mass, fi. As already stated above, we will 
largely concentrate on the /i < case. Together with At > this yields At < 0, which 



is a sufficient condition [32| for the MSSM prediction on BR{b — 57) to be compatible 
with experiment in the presence of a relatively light charged Higgs boson (as the one par- 
ticipating in the top decay under study). In fact, it is known that charged Higgs bosons 
of (9(100) GeV interfere constructively with the SM amplitude and would render a final 
value of BR{b —>■ s'j) exceedingly high. Fortunately, this situation can be remedied in 
the MSSM since the alternative contribution from charginos and stops tends to cancel the 
Higgs contribution provided that Atfi<0. Furthermore, one must also require relatively 
light values for the masses of the lightest representatives of these sparticles, as well as 



high values of tan P |^ ; hence one is led to a set of conditions which fit in with nicely to 
build up a favourable scenario for the decay t —>■ b. 

The evolution of the individual contributions (|83|) , together with the total MSSM 
yield, as a function of yU < 0, is shown in Figs.lla-llb for given values of the other 
parameters. We immediately gather from these figures that the SUSY-QCD correction is 
extremely sensitive to /i. In fact, ^susy-qcd grows rather fast with This is already 
patent at the level of the leading 6mb/mb effect given by eq. (fflD . In all figures where a 
definite /i < is to be chosen, we have taken the moderate value /i = — 150GeV^. In 
this way, for M ~ = 150 GeV, we guarantee that the lightest chargino mass is above 



the LEP 1.5 phenomenological bound: m + > 65 GeV 26 . Concerning the electroweak 
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contribution, we noted above that the component 6mb/mb, eq.(^), actually decreases 
with fi. However, the /i dependences in the full ^susy-ew are more complicated than in 
^susY-QCD and cannot be read off eq.(^). This is evident from Fig. 11a where the total 
(5susY-EW is fairly insensitive to fi; i^mssm, therefore, inherits its marked /^-dependence 
basically from the SUSY-QCD component. As for the sensitivity of the corrections on 
the 5'f/(2)i-gaugino soft SUSY-breaking parameter, M, Fig. 12 conveys immediately that 
it is virtually non-existent. We point out that the choice of M and n is always made in a 
range where the chargino masses are above 65 GeV. 

There is some slight evolution of the corrections with Ab (Fig. 13a), mainly on the 
SUSY-QCD side. We realize that ^susy-qcd is not perfectly symmetric with respect to 
the sign of Ab. Once the sign < is chosen, the correction is larger for negative values 
of Ab than for positive values. We have erred on the conservative side by choosing Ab = 
+300 GeV wherever this parameter is fixed. As far as At is concerned, 5susy-qcd can only 
evolve as a function of that parameter through vertex corrections, which are proportional 
to At cot P (Cf. Appendix A); however, at large tan/? these are very depressed. The 
electroweak correction 5susy-eW) instead, is very much dependent on At; indeed, a typical 
component exhibiting this behaviour is given by eq.(|76|), which is linear in At. The full 
dependence, however, is not linear and is recorded in Fig. 13b. We realize that 5susy-ew 
and 5mssm change sign with At. The shaded vertical band in Fig. 13b is excluded by our 
choice of parameters in Fig. 8. 

Another very crucial parameter to be investigated is the value of m^^ . This is because 
the SUSY-QCD correction hinges a great deal on the value of the sbottom masses, as it 
is plain from eq.(0). As a matter of fact, a too large a value of rrij^^ may upside down 
the leadership of the SUSY-QCD effects. As a phenomenological lower bound for all the 



squark masses we take the absolute LEP 1.5 bound rriq > Q5GeV However, as a 
typical mass value for all squarks other than the stop we use rriq > 150 — 200 GeV{q ^ t). 
From Fig. 14a we see that provided m^^ < 300 Gel^ the SUSY-QCD effects remain domi- 
nant, but they steadily go down the larger is m^^. The electroweak correction Ssusy-ew, 
on the other hand, is quite sustained with increasing mi_^ and there are parameter config- 
urations where for sufficiently heavy sbottoms the supersymmetric electroweak effects are 
larger than the SUSY-QCD effects. However, this is not the most likely situation. The 
behaviour of the branching ratio is plotted in Fig. 14b. 

Obviously, the evolution of the SUSY-QCD corrections with the stop masses is basi- 
cally flat (Fig. 15) since the leading contribution is independent of m^-^. Therefore, it is of 
little help to use the strict lower mass bound m^-^ ~ 65 GeV in our calculation instead of, 
say, the more conservative m^-^ ~ 100 GeV. Nonetheless, if we wish to keep ^mssm > 0, we 
cannot go too far with m^^ , for the electroweak correction is also seen to decrease with m^^ . 
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Indeed, whereas for m^-^ = 65 — 100 GeV one has 5susy-ew ~ 20%, for m^^ ^ 250 GeV one 
finds 5susY-EW ~ 10%. For heavier stop masses, (5mssm becomes zero or slightly negative. 
In this situation, the imprint of SUSY lies in the fact that the total quantum effect is not 



as negative as predicted by standard QCD, eg. (|85D . 

The influence from the sleptons and the other squarks is practically irrelevant as 
it is borne out by Figs.l6a-16b. They enter the correction through oblique (universal) 
quantum corrections. The only exception are the r-sleptons ("staus"), since they are 
involved in the process-dependent (non-oblique) contribution eq. (|i8|) , where the r-lepton 
Yukawa coupling becomes enhanced at large tan/5. For this reason, ^susy-ew in Fig. 16b 
is somewhat larger the smaller is the r-sneutrino mass (assumed to be degenerate with 
the other sneutrinos). In all our calculation we have fixed the common sneutrino mass at 
200 GeV^. 

We have also tested the variation of our results as a function of the external particle 
masses, namely the masses of the top quark, bottom quark and charged Higgs. As for the 
external fermion masses, the corrections themselves are not very sensitive (see Figs. 17a 
and 18a). Among the SUSY corrections, the most sensitive one on rrit (respectively on 
mfe) is 5susY-EW (resp. 5susy-qcd)- The branching ratios also show some dependence 
(Figs. 17b and 18b), especially on rrih. This effect is mainly due to the variation of the 
tree-level partial widths as a function of rrit and mf,. As for the charged Higgs mass, 
Mh+, up to now it has been fixed at Mh+ = 120 GeV. We confirm from Fig. 19a that 
there is nothing special in the chosen value for that parameter since the sensitivity of the 
correction is generally low, except near the uninteresting boundary of the phase space 
where the branching ratio (Fig. 19b) boils down to zero. 

We close our study of the corrections in the yU < case by plotting 6^ as a function of 
tan P (see Fig. 20). By definition, 6r is that part of (5mssm originating from the full process- 
dependent term A,-, eq . (^Sf ) , which stems from our definition of tan/5 on eq.(^). This 
piece of information is relevant enough. In fact, it should be recalled that the quantum 
corrections described in the previous figures are scheme dependent. In particular, they 
rely on our definition of tan/5 given on eq.(^5|). What is not scheme dependent, of course, 
is the predicted value of the width and branching ratio (Figs. 8 and 9b) after including all 
the radiative corrections. Now, from Fig. 20 it is clear that the AT--term is not negligible, 
and so there is a pro cess- dependence in our definition of tan/5, as it was announced in 
Section 3. At first sight, the 5T--effects are not dramatic since they are small as compared 
to 5susY-QCD, but since the latter is cancelled out by standard QCD we end up with 6r 
being of the order (roughly half the size) of the electroweak correction (5susy-ew- 

The main source of process-dependent ^-r-effects lies in the corrections generated by the 
r-mass counterterm, 6mr/mr, and can be easily picked out in the electroweak-eigenstate 
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basis (see Fig. 21) much in the same way as we did for the 6-mass counterterm. There 
are, however, some differences, as can be appraised by comparing the diagrams in Figs. 7 
and 21, where we see that in the latter case the effect derives from diagrams involving 
r-sleptons with gauginos or mixed gaugino-higgsinos. An explicit computation of the 
diagrams (a) + (b) in Fig.21 yields 



,2 



/i M' tan P I{mf^ , ruf^ , M') 



+ -|-^/iM tan/? J(/i,mp^,M), (86) 

where g' = gsw/cw and M',M (Cf. Appendix A) are the soft SUSY-breaking Ma- 
jorana masses associated to the bino B and winos W^, respectively, and the function 
I{mi,m2,m^) is again given by eg. ([75|) . In the formula above we have projected, from 
the bino diagram in Fig. 21a, only the leading piece which is proportional to tan/5. Even 
so, the contribution from the wino-higgsino diagram in Fig. 21b is much larger. Numerical 
evaluation of the sum of the two contributions on eq. (^) indeed shows that it reproduces 
to within few percent the full numerical result (Cf. Fig.20) previously obtained in the 
mass-eigenstate basis, thus confirming that eq.(^) gives the leading contribution. In 
practice, for a typical choice of parameters as in Fig.8, this contribution is approximately 
cancelled out by part of the electroweak supersymmetric corrections associated to the 
original process t b, and one is effectively left with eq. ([76|) as being the main source 

of electroweak supersymmetric quantum effects at high tan j3. 

Finally, the corrections corresponding to the case where fi > are studied in Figs. 22a 
and 22b. The problem with /i > is that, then, the large SUSY-QCD corrections have 
the same (negative) sign as the conventional QCD effects, and as a consequence the total 
MSSM correction can easily blow up above 100%, the branching ratio becoming negative!. 
To avoid this disaster (from the point of view of perturbation theory), we enforce the 
SUSY-QCD correction to be smaller than in the /i < case by assuming an "obesse 
SUSY scenario" characterized by very large values for the sbottom mass (m^^ = 600 GeV) 
and the gluino mass {rrig = 1000 GeV). We also choose At < so that the electroweak 
SUSY correction becomes opposite in sign to the SUSY-QCD correction (a feature that 
also applies in the /i < case, see Fig. 13b) and in this way the total SUSY correction 
is further lessened in absolute value. Incidentally, we remark that the simultaneous sign 
change of both /i and At is also necessary in order to keep Atfi<0; as noted above, this is 
required in order that the MSSM can be compatible with BR{b ^ S7) in the presence of 
a relatively light charged Higgs. In Fig. 23 we bring forward the effect of the new situation 
on the total partial width. In the present instance, the physical signature would be to 
measure a partial width significantly smaller than the one predicted by QCD. Clearly, the 



38 



H > {At yU < 0) scenario is not as appealing as the /i < {At H < 0) one. 

In the end, from the exphcit numerical analysis, we have confirmed our expectations 
that the SUSY-QCD contribution to r{t H^h) is generally dominant. This conclusion 
would not hold only in some (unlikely) cases, e.g. if the gluino is very light and/or the 
lightest bottom squark is "obesse" as compared to lightest top squark, i.e. if the former is 
unusually much heavier than expected. Furthermore, by restricting ourselves to the case 
/i < {At < 0) we confirm that at large tan (3 and for typical values of the parameters 
the total (standard plus supersymmetric) QCD correction largely cancels out, leaving 
a remainder on the SUSY-QCD side (Figs. 8-9). In all circumstances the virtual Higgs 
effects remain comparatively very small. Around tan/? = rrit/mh — 35, one is left with 
basically the electroweak supersymmetric correction, Ssusy-ew: which can be sizeable 
enough to be pinned down by experiment. However, as stated above, there is in general 
a strong remainder, Ssusy~qcd + Sqcd > 0, which grows very fast with tan/5 and it has 
the same sign as Ssusy-ew- In this favourable scenario, the virtual SUSY effects could be 
spectacular. This is true not only because in the relevant window of parameter space the 
SUSY quantum corrections are by themselves rather large, but also because they push 
into the opposite direction than the "expected" standard QCD corrections. As a result, 
the relative deviation between the MSSM prediction and the QCD prediction effectively 
"doubles" the size of the observable effect, a fact which is definitely welcome from the 
experimental point of view. 

From all the previous discussion there is one fact standing out which can be hardly 
overemphasized: If the charged Higgs decay mode of the top quark, t b, does show 

up with a branching ratio of order 10% or above (perhaps even as big as 50%), a fairly 
rich event statistics will be collected at the Tevatron and especially at the LHC e.g. by 
making use of the identification methods described in Section 2. If, in addition, it comes 
out that the dynamics underlying that decay is truly supersymmetric, then the valuable 
quantum signatures that our calculation has unveiled over an ample portion of the MSSM 
parameter space should eventually become manifest and, for sure, we could not miss them. 

At present all the collected event statistics basically relies on our experimental ability 
to recognize the top quark decays originating from standard patterns (angular distribu- 
tion, energy spectrum, jet topology etc.) associated to the usual Drell-Yan production 
mechanism. Notwithstanding, we wish to point out that it should in principle be pos- 
sible to clutch at the supersymmetric virtual corrections associated to the vertex t b 
also through an accurate measurement of the various inclusive top quark and Higgs bo- 
son production cross sections in hadron colliders. As an example, in Fig. 24 we sketch a 
few alternative mechanisms which would generate top quark production patterns heavily 
hinging on the properties of the interaction t b H^-vertex . Thus, while this vertex 
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could be responsible in part for the decay of the top quark once it is produced, it might 
as well be at the root of the production process itself at LHC energies, where it could 
take over from Drell-Yan production |^ . 



We observe that in some of these mechanisms a Higgs boson is produced in associa- 
tion, but in some others (fusion processes) the Higgs boson enters as a virtual particle. 
Now, however different these production processes might be, all of them are sensitive to 
the effective structure of the tb H^-veitex. Similar mechanisms can of course be depicted 
involving the neutral Higgs bosons of the MSSM interacting with ti and bb via enhanced 
Yukawa couplings While it goes beyond the scope of this paper to compute the SUSY 
corrections to the production processes themselves, we have at least faced the detailed 
analysis of a partial decay width which involves one of the relevant production vertices. 
In this way, a definite prediction is made on the properties of a physical observable and, 
moreover, this should suffice both to exhibit the relevance of the SUSY quantum effects 
and to demonstrate the necessity to incorporate these corrections in a future, truly com- 
prehensive, analysis of the cross-sections, namely, an analysis where one would include the 
quantum effects on all the relevant production mechanisms within the framework of the 
MSSM. For this reason we think that in the future a precise measurement of the various 
(single and double) top quark production cross-sections @, ^ should be able to detect 
or to exclude the tb H^-veitex as well as the vertices qq A^{h^, H^) involving the neutral 
Higgs particles of the MSSM and the third generation quarks q = t,b. 

We conclude our discussion with the following remark. Whereas, on the one hand, one 
expects that some top quark partial widths will be determined with an accuracy of 10% 
at the upgraded Tevatron and perhaps better than 10% at LHC ||2^, on the other hand 
we believe that from the point of view of an inclusive model-independent measurement 
of the total top-quark width, Tt, the future e"*" e~ supercollider should be a better suited 
machine |]66| , |67| . For, in an inclusive measurement, all possible non-SM effects will appear 
on top of the corresponding SM effects already computed in the literature |5^. Moreover, 
as shown in Ref.||66|, one hopes to be able to measure the total top-quark width in e"*" e~ 
supercolliders at an unmatched precision of ~ 4% on the basis of a detailed analysis of 
the threshold effects in the cross-section, in particular of the top momentum distribution 
and the resonance contributions to the forward-backward asymmetry in the ti threshold 
region. Under the assumption that F// ~ T^y, and that the SUSY effects on F^ are purely 
virtual effects, it follows that a large SUSY correction of, say 50%, to t ^ b translates 
into a 20% correction to F^. This effect could not scape detection. Thus, the combined 
information from a future e"*" e~ supercollider and from present and medium term hadron 
machines can be extremely useful to pin down the nature of the observed effects. Our 
general conclusion is, therefore, extremely encouraging: In view of the potentially large 
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size and large variety of manifestations, quantum effects on top quark and Higgs boson 
physics could be the clue to the discovery of "virtual" Supersymmetry. 
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Appendix A. Interaction Lagrangian of the MSSM in the 
mass-eigenstate basis 

In this appendix we single out specific pieces of the MSSM Interaction Lagrangian 
involved in the decay rate of the process t b. Although the Lagrangian of the 

MSSM is well-known PI , it is always useful to project explicitly the relevant pieces and 
to cast them in a most suitable form for specific purposes. As a matter of fact, we have 
produced a complete set of Feynman rules for the MSSM using an algebraic computer 
code based on MATHEMATICA ||68|| p^. Here we limit ourselves to quote the Lagrangian 
interactions affecting the process-dependent parts of our decay and omit the interaction 
pieces needed to compute the universal counterterm structures, whose SM parts are well- 



known |^0| and the corresponding SUSY contributions are also available in the literature 
since long time ago ||4^ . All our interactions are expressed in the mass-eigenstate basis. 
Within the context of the MSSM we need two Higgs superfield doublets 



^^We have corrected several mistakes in the subset of rules presented in Ref.|Q involving sparticles 
and Higgses. 
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with weak hypercharges Yi^2 = =Fl- The (neutral components of the) corresponding scalar 
Higgs doublets give mass to the down (up) -like quarks through the VEV < >= Vi 
(< H2 >= V2). This is seen from the structure of the MSSM superpotential 

W = ei, [K H{QW + hMQ'lJ - iJiHlHi] , (A.2) 

where we have singled out only the Yukawa couplings of the third quark-squark generation, 
(t, h) — {t, b), as a generical fermion-sfermion generation of chiral matter superfields Q, U 
and D. Their respective scalar (squark) components are: 



Q= I' , U = i'l , D = r^, (A.3) 



b' 



L 



with weak hypercharges Yq = +1/3, Yu = —4/3 and = +2/3. The primes in ( [A.3[ ) 



denote the fact that q'^ = {q'l^, q'^} are weak-eigenstates, not mass-eigenstates. The ratio 

tan/?=— , (A.4) 

Vl 

is a most relevant parameter throughout our analysis. 

We briefly describe the necessary SUSY formalism entering our computations: 

• The fermionic partners of the weak-eigenstate gauge bosons and Higgs bosons, 
called gauginos, B, W, and higgsinos, H, respectively. From them we construct 
the fermionic mass-eigenstates, the so-called charginos and neutralinos, by, first, 
forming the following three sets of two-component Weyl spinors: 

Tt = {-tW^,Hn, rr = {-tW-,H{} , (A.5) 

Tl = {-zBO,-zWlHlHf} , (A.6) 

which get mixed up when the neutral Higgs fields acquire nonvanishing VEV's and 
then diagonalizing the resulting "ino" mass Lagrangian 

/ M V2MwCp\. 



Cm = - < r+i ( ---w-i^ ) ir- > 



( M' Mzspsw -Mzcpsw\ 

M -MzSf}Cw Mzc/scw 

MzSpSw -MzSpCw -/X 

\-MzCpSw MzCfsCw -/i / 

+ h.c. , (A.7) 



- -<ro 

2 



where we remark the presence of the parameter fi introduced above and of the 
soft SUSY-br caking Majorana masses M and M', usually related as M'/M = 
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(5/3) tan^6'vK, and where cp = cos/3 and sp = sin/3. The corresponding mass- 
eigenstates[3 (charginos and neutrahnos) are the following: 

"^t = ( ) , ^7= C^^~^ = ( ) , (A.8) 



and 

where the matrices N are defined through 

U*MV^ = diag{Mi, M2} , N*M^N^ = diag{M^, . . . , M°} . (A.IO) 

In practice, we have performed the calculation with real matrices U, V and A^, so we 
have been using unphysical mass-eigenstates (associated to non-positively definite 
chargino-neutralino masses). The transition from our unphysical mass-eigenstate 
basis {\E'} = {\E'f,\E'°} into the physical mass-eigenstate basis {x} = {xf^Xa} can 
be done by introducing a set of e parameters as follows: for every chargino-neutralino 
whose mass matrix eigenvalue is Mj, Mq, the proper physical state, is given by 



^ 11 t = ± / A 1 1 \ 




and the physical masses for charginos and neutralinos are m^± = e Mj and m^o = 
eM°, respectively. Needless to say, in this real formalism one is supposed to propa- 
gate the e parameters accordingly in all the relevant couplings, as shown in detail in 
Ref . 1^ . This procedure is entirely equivalent ||6^ to use complex diagonalization 



matrices insuring that physical states are characterized by a set of positive-definite 
mass eigenvalues; and for this reason we have maintained the complex notation in 
all our formulae in Section 4. Whereas for computations with real sparticles the 
distinction matters |l23[| , for virtual sparticles the e parameters cancel out, and so 



one could use either basis {\E'} or {x} without the inclusion of the e coefficients. We 
have stressed here the differences between the two bases just to make clear what are 
the physical chargino-neutralino states, when they are referred to in the text. 

Among the gauginos we also have the strongly interacting gluinos, g"^ {r = 1, . . . ,8), 
which are the fermionic partners of the gluons. 

As for the scalar partners of quarks and leptons, they are called squarks, q, and slep- 
tons, I, respectively. Again we will use the third quark-squark generation (t, b) — 
{t, b) as a generical fermion-sfermion generation. The squark mass-eigenstates. 



^^We use the following notation: first Latin indices a,b,...=l,2 are reserved for sfermions, middle Latin 
indices i,j,...=l,2 for charginos, and first Greek indices a, /3, ... = 1, . . . , 4 for neutralinos. 
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Qa = {(11,(12}, if we neglect intergenerational mixing, are obtained from the weak- 
eigenstate ones q'^ = {q'l = q^, q'2 = qn}, through 

b 

= -^^^/O iq = t,b). (A.12) 

\sm On cos Oa J ' ' 



The rotation matrices in ([A .121) diagonahze the corresponding stop and sbottom 
mass matrices: 

fMl+ml + cos 2/3(Tf - Q,.2^)M| m,Ml^ 
' V rrigMlj, Ml + + cos 2/5 s^, M| 

(A.13) 

Ri'i)iMlR^'^ = dtag{ml,ml} (m^, > rn^-J , (A. 14) 

with T^^ the third component of weak isospin, Q the electric charge, and Mg^ ^ the 
soft SUSY-breaking squark masses [|[. (By SU (2) i-gauge invariance, we must have 
Mf^ = M^^, whereas M^^, M^^ are in general independent parameters.) The mixing 
angle on eg. ( [A. 121) is given by 



t^nZU, - MI-MI + cos2/3(Tr - 2g,.^)M| ' ^""'''^ 

where 

M|^^ = A - cot /3 , Mlji = Ah- fitanp , (A.16) 

are, respectively, the stop and sbottom off-diagonal mixing terms on eq. (|A.13| ). 
Furthermore, /i is the SUSY Higgs mass parameter in the superpotential, and At^ 
are the trilinear soft SUSY-breaking parameters. We shall assume that \At^b\ < 3Mg, 
where Mq is the average soft SUSY-breaking mass appearing in the mass matrix 
(|A.13D; this relation roughly corresponds to the necessary, though not sufficient. 



condition for the absence of colour-breaking minima . 

The charged slepton mass-eigenstates can be obtained in a similar way after straight- 
forward substitutions in the mass matrices, with the only proviso that there is no upt, 
so that ul is itself the sneutrino mass-eigenstate, hence R^^J = unless a = b = 1 
where R^^i = 1. 

Next let us describe the relevant pieces of the MSSM interaction Lagrangian involving 
the fields defined above. 

• fermion— sfermion— (chargino or neutralino) 
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After translating the allowed quark-squark-higgsino/gaugino interactions into the 
mass-eigenstate basis, one finds 



a=l,2 i=l,2 

+ 4 E E {-i:^l{AtaPL + A^l^Pj,)t + b:n{A^^aPL + A^'^^^^^ 

a=l,2a=l,...,4 

+ h.c. (A.17) 

wViprp /l'-*'' A^''^ A^!^ A^!'^ arp 
Wneie ■ri±ai, ^±ai) ^ztaa; ^±aa die 

A-ai = ~\Rla , 

= (Ar*2 + FLtan^^Ar*J + y2At4?iV*3, 

A^l\a = V^At-Ria*iVo3 — Y^tan6'vK-R2a*A'al , ( A T Q\ 

4(6) _ pW*T/* X p(6)*T/* ^^-^^^ 

^+ai — -Kla Vil - M^2a ^i2 1 

^-ai — ~^t-ttia Ui2 , 

A^'L = -y2Abi?ff*Ar,4 + F^tan^^4f AT^i. 

with Yl and Y^''' the weak hypercharges of the left-handed SU (2)^ doublet and right- 
handed singlet fermion, and and Xf, are - Cf. eq.(0) - the potentially significant 
Yukawa couplings normalized to the SU (2) ^ gauge coupling constant g. 

quark— squark—gluino 

^9,^ = -^€,k 9r {Ria*PL " R^£* Pr) Qi + h-c (A.19) 

where A'' are the Gell-Mann matrices. This is just the SUSY-QCD Lagrangian 
written in the squark mass-eigenstate basis. 



• squark— squark— Higgs 



where we have introduced the coupling matrix 

^ah — J^cb ^da Hcd 

_ / tan /? + cot /3 — sin 2/? m;, (/i + tan/3) \ (p^2\) 
^""^ ~ \ mt{n + At cot (3) mtrrib (tan /? cot /?) / ' ^ ' ' 

chargino—neutralino— charged Higgs 

£^±^^^0 = -gH-^l (cos/3 g^,Pi + sin/? Q^^Pr) ^+ + h.c. (A.22) 

with 

/ Q^. = U;,Nl, + ^ {Nl, + tan e^N^,) U*, 

\ nR — ■[/ AT i fAT I ^-„v./3___Ar [J\.^6) 



Qai = VaN^^-j^{N^2 + tanewN^i)V, 



i2 ■ 
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• Gauge interactions 

In our calculation we only need the sparticle interactions with the W^: 
-queirks 



-squeirks 



^w±,, = ^ trPib W+ + h.c. (A.24) 



jO-w^m = ^^Ril*RwKt: k + h.c. (A.25) 



-charginos and neutralinos 



Appendix B. The Counterterm 5Zhw 

We discuss the structure of 5Zhw both in the unitary gauge and in the Feynman 
gauge. In order to renormalize the mixed self-energy a mixing counterterm 

must be introduced in the wave function of W^: 

Wi^{zffhv-±,'-^3,H-. (B.1) 

B.l Uniteiry gauge 

If one chooses to work in the unitary gauge, then the renormalization is straightfor- 
ward: 

Lug = -\f^.uF^' + M^W^W'^ ^ L,t = tMw6ZHw{W^d^H+ - W^d'^R-) . (B.2) 
In this gauge the corresponding renormalized 2-point Green function reads (Fig. 25a): 



w 



:b.3) 



Thus, a renormalized self-energy can be defined as follows: 

±Hw{k^) = ^Hw{k^) - M^SZhw (B.4) 

Now we must impose a renormalization condition on 'E^^{k'^); and we choose it in a way 
that the physical Higgs does not mix with the physical W^: 
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B.2 Feynman gauge 

As we carry out our calculation in the Feynman gauge, we would also like to perform 
the renormalization of the Higgs sector in that gauge. The Lagrangian is sketched as 
follows 0]: 

L = Lc + Lgf + Lpp. (B.6) 
where Lgf stands for the gauge-fixing term in that gauge, 



-F+F~ + ... {F^ = d^'W;^ T iMw G^) (B.7) 



and Lpp = fj"" {^dF"" / 89'^^ rf is the Faddeev-Popov ghost Lagrangian constructed from 
FP and anti-FP Grassmann scalar fields 77, fj. Since we are interested in the charged 
gauge-Higgs {W^ — H^) and charged Goldstone-Higgs {G^ — H^) mixing terms in that 
gauge, we have singled out just the relevant term on eq. (|B.7|) . 

The relationship with the original weak-eigenstate components on eq.(Al) is the fol- 
lowing: 

As is well-known, although the classical Lagrangian, Lc, also contains a nonvanishing 
mixing among the weak gauge boson fields, W^, and the Goldstone boson fields, G^, 
namely 

= iMw W-d^'G^ + h.c. , (B.9) 

the latter is cancelled (in the action) by a piece contained in Lgf- Now, after substituting 
the renormalization transformation for the Higgs doublets, eq.(^), on the Higgs boson 
kinetic term with SU{2)l x f/(l)y gauge covariant derivative, one easily projects out the 
following relevant counterterms 

5C = dZnid^H+d^'H- + 5ZG±d^,G+d^'G- 



+ 6Zhg [d^H-d''G^ + h.c.j + 6Zhw [iMw W-d''H+ + h.c.j + ... (B.IO) 

where 



and 



SZhg = spcp {5Zh2 - 5Zh^ 



{oZh2 - oZhJ + 



2' tan/? 
The renormalization transformation for the VEV's of the Higgs potential 



(B.12) 



V, ^ Z]l^{v, + 5v,) = h + ^ + i 5ZhA V, , (B.13) 
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implies that the counterterm to tan/3 is related to the fundamental counterterms in the 
Higgs potential by 

/) tan /3 /)7)r, /)7)i 1 

:b.i4) 



(5 tan/3 5v2 5vi 1 
tan/5 f2 Vi 2 



If one imposes the usual on-shell renormalization conditions for the A^-boson, one has 

tan (3 + cot (3 



6Zh2 - 5Zh^ 



^az{Mao, 



(B.15) 



In addition, there is another mixing term between and originating from the mass 
matrix of the Higgs sector [^|. This one loop mixture is contained in: 



l\/fb2 *o 



H 

G 



_6 



:b.i6) 



where we have attached a superscript b to bare quantities, and tj are the tadpole coun- 



terterms 



to 



:B.17) 



— sin {l3 — a) t^o + cos (/3 — a) t^o 
sin (/? — a) tfiO + cos (/? — a) tuo . 

We are now ready to find an expression for the mixed 2-point Green functions (Figs. 25a 
and 25b). As the mixing Lagrangian between and on eq. (|B.2| ) is the same as on 
eq. (|B.10|) , the 2-point Green function will have the same expression (|B.3|) but with the 
propagator in the Feynman gauge, namely: 



k' 



.-iT^Hwik'^) 



+ ik" M? 



6Z 



HW 



W 



W 



w 



k^-M^ 



k^-Mjj^ L 



ki, 



-itHw{k^) 



Ml 



w 



k^-M^ 



A 



HW 



:b.i8) 



Next we impose Si^iy(M^±) = as before, leading to the same formal expression for 
SZhw as in eq.( p.5| ). However, as a new ingredient, we now have the mixed if^ — 
self-energy: 



k^ - Ml± 



This allows us to define renormalized self-energies. 



and 



:b.i9) 



-'HG 



-'HQ 



kHZHG + 



(B.20) 



The mixed self-energies T,Hw{k'^) and T^Haik'^) obey the following Slavnov- Taylor identity: 



(B.21) 
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This identity is derived from a BRS transformation involving the Green function con- 
structed with an anti-FP field and the charged Higgs field: < 0\6brs {v^ H^)\0 >= 0. 



Following the standard procedure [|7T| one immediately gets 



< 0|F+ H+\0 >=< Old^W^ 
which in momentum space reads 



H 



iMwG+H+\0 >= 0, 



HW 



with 



Mw A^^ = 



k^ 



:B.22) 



(B.23) 



(B.24) 



Clearly, eq. (p.23| ) implies eq. ( |B.21| ) . The latter identity guarantees that the contribution 
to our decay t —>■ b from the two diagrams in Figs. 25a and 25b vanishes since no 
mixing is generated among the physical boson and the renormalized fields and 
W^: T,hg{M]j±) = T,Hw{Mfj±) = 0. There is of course another Slavnov- Taylor identity, 
derived in a similar manner, which insures that the renormalized mixing between and 
also vanishes. Thus we have proven that the expression for 6Zhw is formally the 
same in both unitary and Feynman gauges, but that in the latter gauge one must take 
into account the additional renormalization of the mixed self-energy T^hg- 



Appendix C. Vertex functions 

In this appendix we briefly collect, for notational convenience, the basic vertex func- 
tions frequently referred to in the text. Dimensional regularization is used throughout. 
The given formulas are exact for arbitrary internal masses and external on-shell momenta. 
Most of them are an adaptation to the g^^i, = {H } metric of the standard formulae 



of ref.|72]. The basic one, two and three-point scalar functions are: 



Aoim) = I rq ^ (C.l) 

G,{p, k, m„ m„ ma) = / d^g _ + [(, + ^ + _ ' (C-3) 

using the integration measure 

rq = ^j,(^-^)-^. (C.4) 
The two and three-point tensor functions needed for our calculation are the following 
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[Co,Cf„C^u]{p,k,mi,m2,m3) = 

By Lorentz covariance, they can be decomposed in terms of the above basic scalar func- 
tions and the external momenta: 

Bf,{p,mi,m2) = Pi^Bi{p,mi,m2). 

B^^{p,mi,m2) = Pf^PuB2i{p,mi,m2) + gf,uB22{p,mi,m2). 

Co{p,k,mi,m2,ms) = Bo{k,m2,m:i) + mjCoip, k,mi,m2,ms). 

C^^{p, k, nil, m2, mg) = PnPuC2i + k^k^C22 + {p,xku + k^p^)C23 + 9,xuC2i , (C.7) 

where we have defined the Lorentz invariant functions: 



Bi{p,mi,m2) = ^[Aolmi) - Ao{m2) - fiBo{p,mi,m2)], (C.8) 
Zp 



B2i{p,mi,m2) = 77^7 -^[{n - 2)Ao{m2) - 2mjBo{p,mi,m2) 

Ip-^yn — 1) 

- nfiBi{p,mi,m2)], (C.9) 
B22{p,mi,m2) = -^^^j—^[AQ{m2) + 2mlBo{p,mi,m2) + fiBi{p,rni,rn2)],{C.10) 

( ^^^\^y( Bo{p + k, mi, ms) - Bo{k, 1712,1713) - fiCo \ (C 11) 

\ Ci2 ) \ Bo{p,mi,m2) - Bo{p + k,mi,m3) - f2Co y ' 

( = Y ( + ^' "^1' "^3) + Bo{k, ma, mg) - fiCu - 2C24 \ -^2) 

\ C23 J \ Bi{p,mi,m2) - Bi{p + k,mi,m3) - f2Cii )' 

+ k,mi, 1713) - /2C12 - 2C24]} , 

C24 = — ^— [5o(A;, m2, ms) + 2m?Co + /idi + /2C12] , (C.14) 
2[n — 2) 

the factors /i_2 and the matrix y, 

r 2,2 2 

/I = P +1711 -^2, 

/2 = /c^ + 2pk + ml — ml , 

^ " 2[p2p - (pky] ( -pA; ) ■ (C.15) 

The UV divergences for n — > 4 can be parametrized as 

e = n — 4, 

A = ^ + ^^-ln(47r), (C.16) 
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being 7^ the Euler constant. In the end one is left with the evaluation of the scalar 
one-loop functions: 



Ao{m) 
Bo{p,mi,m2) 



IGvr^ / ' ,,2 



m^(A - 1 + In — ) , 

2 

A + In ^ - 2 + ln[(xi - l){x2 - 1)] 



, 1 ^1^1 ^2 
-l-Xi m + X2 m 



xi — 1 

i \l 1 



X2 - 1 



with 



/ \ TTZi ^^2 I \ 1 /2 / 2 2 2\ 

2^1,2 = a:;i,2(p,mi,m2) = -H — ± —\ > [p ,m^,m^), 

X(x,y,z) = [x - (^/y - y/zy][x - (y/y + y/zf] , 



(C.17) 

(C.18) 
(C.19) 

(C.20) 



and where S is a bookkeeping device for the following alternate sum of twelve (complex) 
Spence functions: 



E = Sp[ 
—Sp 
Sp( 



yi 



( 

\y2- zf 

ys - zf) 



-Spy 

-Sp[ 



V\ - 




yi - 


z\) 


( y2 


- 1 


\y2 - 


~ii 
~ Zi 


ys - 


-1 ^ 


ys- 


-ill 

^1 y 



+ Sp 



yi 



.y^ - ^2, 
( y2 



Sp 



\y2 - ^2 , 

ys 

M - zf\ 



Sp 
+ Sp 
— Sp 



yi - zi) 

' y2-i \ 
yy2-z^) 

' 1/3 - 1 ' 
,y3 - zf ^ 



(C.21) 



The Spence function is defined as 

Sp{z) 

and we have set, on one hand: 



1 ln(l - zt) 



dt, 



(C.22) 



^1,2 



"1,2 



Zl,2 — 



a^i,2(p, m2,mi) , 
Xi,2{k,m3,m2) ; 



and on the other: 

yi = yo + C , y2 

where 



yo 



yo 



yo = -7; 



2 pk + ' 



g — —k'^ + ml — ml, h — —p^ — 2pk — ml + ml, 



(C.23) 



(C.24) 



(C.25) 
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and ^ is a root (always real for external on-shell momenta) of 

+ 2pki + /c^ = . 



(C.26) 



Derivatives of some 2-point functions are also needed in the calculation of selfenergies, 
and we use the following notation: 

d_ 

Qp2 



^B^{p, nil, 1112) = Bl{p,mi,m2). 



(C.27) 



We can obtain all the derivatives from the basic B'r,: 



5o(p,mi,m2) 



-i \ (]_ 1 
167rV \p2 + Xi/^{p^,mlm^2) 
Xi — 1 



Xi{xi — 1) In 



Xi 



X2[X2 - 1) In 



X2- I 
X2 



, (C.28) 



which has a threshold for \p\ — ■mi + 777.2 and a pseudo-threshold for \p\ — |777i — 7772 1. 
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Figure Captions 

Fig.l The lowest-order Feynman diagram for the charged Higgs decay of the top 
quark. 

Fig. 2 Feynman diagrams, up to one-loop order, for the electroweak SUSY vertex 
corrections to the decay process t — > b. Each one-loop diagram is summed 
over all possible values of the mass-eigenstate charginos {"^f ;i = 1,2), neutralinos 
(\E'° ; a = 1, 4), stop and sbottom squarks {ba, ib]a,b = 1,2). 

Fig. 3 Feynman diagrams, up to one-loop order, for the Higgs and Goldstone boson 
vertex corrections to the decay process t — ^ b. 

Fig. 4 Electroweak self-energy corrections to the top and bottom quark external 
lines from the various supersymmetric particles, Higgs and Goldstone bosons. 

Fig. 5 Corrections to the charged Higgs self-energy from the various supersymmetric 
particles and matter fermions. Only the third quark-squark generation is illustrated. 

Fig. 6 Corrections to the mixed W'^ — H^ self-energy from the various supersymme- 
tric particles and matter fermions. Only the third quark-squark generation is illus- 
trated. 

Fig. 7 (a) Leading SUSY-QCD contributions to 6mb/mb in the electroweak-eigenstate 
basis; (b) Leading supersymmetric Yukawa coupling contributions to Smb/rrib in the 
electroweak-eigenstate basis. 

Fig. 8 The total partial width, TMssMit b), including all MSSM effects, versus 

tan/?, as compared to the tree-level width and the QCD-corrected width. Also 
plotted is the tree-level partial width of the standard top quark decay, t b. 
The masses of the top and bottom quarks are rrit = 175 GeV and = 5GeV, 
respectively, and the rest of the inputs are explicitly given. We remark that At = 
Af, = ... = A is a common value of the trilinear coupling for all squark and slepton 
generations. Unless explicitly stated otherwise, the inputs staying at fixed values in 
the remaining figures are common to the values stated here. 

Fig. 9 (a) The relative corrections 6, eq.(^), as a function of tan/?. Shown are 
the SUSY-EW, standard EW (i.e. non-supersymmetric electroweak), SUSY-QCD, 
standard QCD, and total MSSM contribution, eq. (|83|) ; (b) The branching ratio 
(|79|), as a function of tan/3; separately shown are the values of this observable 
after including standard QCD corrections, full MSSM corrections, and the tree-level 
value. 
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Fig. 10 (a) The relative corrections 6, eg. ([77|) , as a function of the gluino mass, 
rrig, for the SUSY-EW, standard EW, SUSY-QCD, standard QCD contributions, 
and total MSSM contribution, (b) As in (a), but for the tree- level and full MSSM- 
corrected branching ratios (^). We have set tan/5 = mt/rrih = 35. This value is 
maintained in all figures where tan /3 is fixed. 

Fig. 11 (a) The relative corrections S, eg. ([77|) , as a function of the supersymmetric 
Higgs mixing parameter, /i, for the various contributions as in Fig.9a; (b) As in (a), 
but for the same branching ratios as in Fig. 10b. 



Fig. 12 Dependence of 6, eg . ([TTD , on the 5'f/(2)L-gaugino soft SUSY-breaking mass, 
M, assuming that the U{1)y gaugino mass, M', is related to M through M' /M = 
I tan^^vK- The same individual and total contributions as in Fig.Qa are shown. 

Fig. 13 (a) The relative corrections 6, eg . ([77|) , as a function of the trilinear soft 
SUSY-breaking parameter Af, in the bottom sector. The other trilinear couplings 
are kept as in Fig. 8; (b) As in (a), but for the trilinear soft SUSY-breaking parameter 
At in the top sector. Shown are the same individual and total contributions as in 
Fig.9a. 



Fig. 14 (a) The relative corrections S, eg.(|77D, as a function of the lightest sbottom 



mass, m^^, for the various contributions as in Fig. 9a; (b) As in (a), but for the same 
branching ratios as in Fig. 10b. 

Fig. 15 The relative corrections S, eg. ([77|) , as a function of the lightest stop mass, 
m^^, for the various contributions as in Fig.9a. 



Fig. 16 (a) The relative corrections 6, eg . (|77|) , as a function of the up-sguark masses 



niu^ . The c-sguarks are assumed to be degenerate with the up-sguarks; 
(b) (5 as a function of the sneutrino masses, assumed to be degenerate. 

Fig. 17 (a) The relative corrections 6, eg. (|77D , as a function of the top guark mass 
within about 2 a of the present experimental range at the Tevatron. (b) As in (a), 
but for the same branching ratios as in Fig. 10b. 

Fig. 18 (a) The relative correction S, eg.(^), as a function of the bottom guark 
mass, (b) As in (a), but for the same branching ratios as in Fig. 10b. 

Fig. 19 (a) The relative corrections 5, eg.(^), as a function of the charged Higgs 
mass; (b) As in (a), but for the same branching ratios as in Fig. 10b. 



59 



Fig. 20 The supersymmetric (5tsusy-ew) and non-supersymmetric (^tew) elec- 
troweak contributions to 6, eq.([77D, from the pro cess- dependent term A^, eq. (^8D , 
as a function of tan /3. 

Fig. 21 Leading supersymmetric electroweak contributions to 8171^-/771^ in the electro- 
weak-eigenstate basis. 

Fig. 22 (a) The relative corrections 5, eq.(|77D, as a function of the hghtest sbottom 
quark mass, for positive /i = +150 GeT^, negative At and given values of the other 
parameters. In this case, huge values of m^^ are needed in order to damp the 
absolute value of the total correction down to 100%; (b) As in (a), but for the same 
branching ratios as in Fig. 10b. 

Fig. 23 The total partial width, V{t H^b), including all MSSM effects, versus 
tan/3, for the same inputs as in Fig. 22a, as compared to the tree-level width and 
the QCD-corrected width. A typical value of m^^ within the range used in Fig. 22a 
is selected. Also plotted is the tree-level partial width of the standard top quark 
decay t — > b. 

Fig. 24 Typical diagrams for top quark and charged Higgs production in hadron 
colliders involving the relevant t b H^-vertex. 

Fig. 25 The mixed blobs — and — at any order of perturbation 
theory. 
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